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ABSTRACT 


This  thesis  presents  a  scheme  for  correcting  certain 
fundamental  errors  in  the  perturbed  stationary  states 
theory  of  slow  atomic  collisions.  The  corrections  are 
associated  with  electron  translation  factors,  of  which  the 
pss  theory  gives  no  account.  In  the  corrected  theory, 
switching  functions  are  used  to  describe  electron  trans¬ 
lation  in  a  molecular  state.  One  can  use  a  single  function 
for  all  electronic  states,  or  each  state  may  have  its  own, 
characteristic  switching  function.  However,  the  functions 
are  not  fully  specified  by  the  theory,  and  the  choice  of 
suitable  ones  is  an  important  practical  problem.  This 
thesis  examines  one  method  of  choosing  them,  based  on  the 
fact  that  the  corrected  nonadiabatic  couplings  for  ionising 
transitions  are  very  sensitive  to  the  choice  of  switching 
function.  For  each  bound  electronic  state,  it  is  possible 
to  choose  an  "optimum"  switching  function,  such  that  the 
corrected  couplings  are  very  much  smaller,  and  of  shorter 
range,  than  those  predicted  by  pss  theory.  These  results 
are  compared  with  those  obtained  using  simpler  choices. 

Ionisation  matrix  elements  are  presented  for  the 

system  of  one  electron  and  two  nuclei,  notably  for  and 

++ 

HeH  ,  and  the  parameters  specifying  the  switching  func¬ 
tions  for  the  lowest  bound  states  are  given.  The 
implications  of  these  results  for  nonadiabatic  couplings 
in  general  are  also  discussed. 
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CHAPTER  I 


INTRODUCTION 

A.  Aims  and  Objectives 

This  thesis  presents  a  scheme  for  correcting  certain 

fundamental  errors  in  the  well-known  and  widely  used 

perturbed  stationary  states  (pss)  theory  of  slow  atomic 

collisions.  The  corrections  in  question  are  associated 

with  electron  translation  factors ,  of  which  pss  theory 

gives  no  account.  The  corrected  theory  has  been  developed 
12  3—6 

by  Thorson  '  and  others  and  has  recently  been  formal- 

7 

ised  by  Thorson  and  Delos,  who  use  switching  functions  to 

describe  electron  translation  in  a  molecular  state.  These 

functions  are  not  fully  specified  by  the  theory,  and  the 

choice  of  suitable  ones  presents  an  important  practical 

problem.  This  thesis  explores  one  method  of  choosing  them, 

based  on  the  fact  that  nonadiabatic  couplings  (corrected 

for  electron  translation  effects)  are  extremely  sensitive 

to  the  choice  of  switching  function,  especially  for 

12 

ionising  transitions.  '  I  show  that  it  is  possible  to 
choose  an  "optimum"  switching  function,  characteristic  of 
each  bound  electronic  state,  such  that  the  corrected  coup¬ 
ling  matrix  elements  are  very  much  smaller,  and  of  shorter 


1 


. 


.. 


range,  than  those  predicted  by  pss  theory.  I  also  compare 

3 

these  results  with  those  obtained  using  simpler  choices 
for  the  switching  function.  For  truly  "molecular"  states, 
the  "optimum"  switching  functions  lead  to  significantly 
better  results,  while  for  purely  "atomic"  states,  the 
method  selects  the  simpler  functions  quite  naturally. 

Detailed  ionisation  matrix  elements  are  presented 
for  the  system  of  one  electron  and  two  nuclei,  especially 
for  and  HeH++,  and  the  parameters  specifying  the 

switching  functions  for  the  lowest  electronic  bound  states 
are  given.  The  implications  of  these  results  for  nonadia- 
batic  couplings  in  general  are  also  discussed. 

In  the  rest  of  the  Introduction,  I  shall  explain  the 
meaning  of  the  terms  used  above,  review  the  history  of  the 
problem,  and  provide  signposts  to  guide  the  reader  through 
the  detail,  which  would  otherwise  obscure  the  main  ideas 
of  the  dissertation. 

B.  Background  ^kterial 

1.  Basics 

The  problem  of  interest  to  us  occurs  in  all  atomic 
and  molecular  collisions,  but  its  essential  features  can 
be  studied  most  easily  in  the  model  systems  considered 
here,  composed  of  two  heavy  particles,  masses  M  ,  Mg,  and 
charges  Z^e ,  ZBe,  and  an  electron,  mass  mQ .  After  separa¬ 
ting  out  the  centre  of  mass  motion,  we  require,  in  prin¬ 
ciple,  the  solution  to  the  time-independent  SchrOdinger 


. 


:  V.  ■ 
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equation , 


H(r,R)  ¥(r,R)  =  E  y(r,R) 


(1-1) 


subject  to  the  appropriate  scattering  asymptotic  boundary 
conditions.  Molecular  co-ordinates,  R,  r,  are  used:  R  is 
the  vector  from  A  to  B,  and  r  is  the  position  vector  of 
the  electron  from  the  centre  of  mass  of  the  nuclei.  For 
the  moment,  we  express  r  in  terms  of  reference  axes  fixed 
in  space,  but,  later  on,  it  will  be  more  convenient  to  use 
a  frame  rotating  with  the  molecular  axis  R.  The  reader 
should  pay  careful  attention  to  the  subject  of  co¬ 
ordinates,  because  it  is  closely  related  to  the  central 
problem  of  the  thesis:  in  an  atomic  collision,  the  system 
separates  asymptotically  into  two  atomic  subsystems,  for 
which  atomic,  rather  than  molecular,  co-ordinates  are  more 
appropriate;  and  the  transformation  of  wavef unctions  from 
one  co-ordinate  system  to  the  other  is  not  quite  as 
trivial  as  it  seems . 


In  molecular  co-ordinates,  the  Hamiltonian,  H(r,R) 


is 


(fi2/2y)  (VR2)  -*  +  he(?;R)  , 


(I-2a) 


v  =  mamb/(W  ; 


(I-2b) 


* 


*  •  ci  *  1 ^  ^  ■  il  I  1  ’■ 


he(r;R)  is  the  molecular  electronic  Hamiltonian, 

he(?;R)  =  -  (fi2/2m)  (^.2)g  +  V(r;R)  ,  (I-3a) 

m  =  mQ  (M^+Mg)/ (M^Mg+niQ)  ,  (I-3b) 

and  V(r;R)  is  the  electrostatic  potential  energy, 

V(r ; R)  =  (ZAZBe2/R)-(ZAe2/rA)-(ZBe2/rB)  .  (1-4) 

For  a  wide  range  of  physical  conditions,  the 
classical  trajectory  approximation , ^ ^  in  which  only  the 
electron  motion  is  treated  quantum  mechanically,  provides 
an  alternative  to  solving  equation  (1-1) .  The  nuclei  are 
assumed  to  follow  some  specified  classical  trajectory  R(t) , 
and  the  electronic  wavefunction  obeys  the  time-dependent 
Schrttdinger  equation, 

h  (r;R(t)  }  T(r,t)  =  ill  {  3T  (r ,  t) /3t  }+.  (1-5) 

We  shall  make  use  of  both  descriptions  of  a  slow  collision 
in  this  thesis;  for  nearly  all  problems  of  interest  to  us, 
the  two  approaches  give  the  same  final  results. 

In  general,  equations  (1-1)  or  (1-5)  are  solved  by 
expanding  the  wavefunction  in  terms  of  some  appropriate 
set  of  electronic  states,  |n(R)>,  projecting  out  the  j'th 


' 
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component  with  <j (R) | ,  and  solving  the  resulting  coupled 
equations.  In  principle,  any  complete  set  of  functions 
will  do,  but,  in  practice,  the  set  is  always  truncated, 
and  the  idea  is  to  choose  a  set  which  approximates  the 
true  behaviour  as  closely  as  possible. 

If  the  nuclear  speed  is  much  less  than  the  electron 
speed,  we  may  expect  that  the  electron  system  will  adjust 
smoothly  to  changes  in  the  nuclear  position;  this  is  the 
fundamental  assumption  of  the  Born-Oppenheimer  approx¬ 
imation  in  molecules,  and  the  same  idea  can  be  used  in 
slow  collision  theory.  At  each  nuclear  configuration,  the 
set  of  Born-Oppenheimer  (fixed-nuclei),  adiabatic  states 
may  be  used  as  the  expansion  basis.  These  states  are 
eigenfunctions  of  the  electronic  Hamiltonian, 

he(r;R)  <|>n(r;R)  =  eR(R)  (j>n(r;R);  (1-6) 

they  depend  parametrically  on  R,  and  coupling  between  them 
occurs  via  this  dependence.  If  the  system  evolves  nearly 
adiabatically ,  coupling  is  unimportant,  and  elastic 
scattering  on  the  potentials  en(R)  dominates.  Conditions 
for  the  validity  of  adiabatic  approximations  are  usually 
discussed  using  the  time-dependent  SchrOdinger  equation; 
it  can  be  shown  that  the  condition  for  adiabatic  evolution 
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|v.<j|^  h  |  n>  |  <<  |  e  (R(t)  }-e  .  (R(t)  }  |/fi,  (1-7) 

i\  “  n  ) 

where  v  =  dR/dt  is  the  nuclear  velocity. 

We  can  now  define  the  term  "slow  collision"  more 
precisely:  a  collision  is  said  to  be  slow  if  the  adiabatic 
criterion  is  satisfied  for  "principal  quantum  transitions" 
--transitions  for  which  As  >  1.0  eV.  For  -  Mg  equal  to 
the  proton  mass,  collision  energies  below  about  one  keV 
are  slow.  However,  note  that,  even  if  v  is  very  small, 
the  adiabatic  criterion  will  be  violated  in  the  case  of 
coupled,  nearly  degenerate  states.  This  means  that  some 
states  of  an  atomic  collision  system  can  never  behave 
adiabatically ,  i.e.,  the  continuum  electronic  states  and 
the  closely  spaced,  Rydberg  levels;  even  for  the  more 
tightly  bound  states,  strong  nonadiabatic  coupling  can 
occur  if  a  local  near-degeneracy  (Rq)  -  ^(Rq)  of  two 
potential  curves  is  encountered. 

Slow  collision  processes  may  thus  be  divided  into 
two  broad  categories : 

a.  Direct  impact  processes ,  for  which  [e^-e^l  is  large, 
and  the  nuclear  momentum  change  accompanying  a  trans¬ 
ition  is  large;  such  transitions  are  improbable,  and 
the  system  evolves  nearly  adiabatically  (e.g.,  impact 
ionisation  of  low-lying,  bound-state  electrons) 

b.  Degeneracy  mediated  processes ,  in  which  the  adiabatic 
criterion  is  locally  violated,  due  to  near-degeneracy 


. 


■ 


II 


■ 


of  two  (or  more)  levels,  and  the  transition  probabil¬ 
ities  are  large,  even  at  low  collision  velocities 
(e.g.,  avoided  crossings  of  adiabatic  potential  curves) 


7 


An  adequate  theory  of  slow  collisions  should  be  able  to 
describe  both  categories. 

2.  PSS  Theory 

The  perturbed  stationary  states  (pss)  method^  takes 
a  truncated  set  of  Born-Oppenheimer  electronic  wave- 
functions  as  the  expansion  basis;  thus,  for  equation  (1-1) , 

N 

¥(?,R)  =  £  xn  (R)  | n  (R)  >  ,  (I-8a) 


|  n  (R)  >  =  <(>n  (r  ;R)  . 


( I  -  8b ) 


One  then  obtains  coupled  equations  for  the  functions  Xr(r) 

.4'. 

as  components  of  an  abstract  vector  _^(R)  : 


( l/2y)  {-ifi  VR  +  P  (R)  }2  +  £(R) 


=  E  )^(R)  . 


(1-9) 


The  matrix,  £(R) ,  is  diagonal,  and  its  nn  element  is  the 
potential  curve,  en(R);  coupling  between  states  j  and  n 
arises  solely  from  the  nonadiabatic  coupling  matrix,  P (R) , 

Pjn(R)  =  -if>  <*jl 


(1-10) 


'  '  ' 


' 


- 
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It  can  be  shown  that  P(R)  is  Hermitian  and  that  its  diag¬ 
onal  elements  vanish.  To  derive  equation  (1-9),  I  have 
used  the  relation 


-ifi  (Vp(R)  +  =  <<(>,  |-fi2(Vp2)i|<t>  >.  (1-11) 

i\  —  —  —  j  ii  j  i\  it  n 


Note  that  the  gradient  in  (1-10)  is  computed  holding  r 
fixed  in  the  non-rotating  frame,  which  accounts  for  the 
existence  of  angular,  as  well  as  radial,  components  in  the 
vector,  P(R)  .  In  actual  calculations,  we  use  a  frame 
rotating  with  the  molecular  axis  to  describe  the  electron 
system;  in  this  frame,  the  Hamiltonian,  hg ,  and  the  basis 


states,  <J>  ,  depend  only  on  R  =  |R[  ,  and  the  angular  parts 

of  P.  arise  directly  from  the  transformation  from  r  (nrf) 
Jn 

to  r'  (rf) .  Physically,  these  terms  represent  the 
Coriolis  interaction. 


In  the  alternative,  classical  trajectory  description, 
the  analogous  pss  expansion  is^ 


T(r,t)  =  £  a  (t)  cf>  (r;R(t)}  exp{-i/te  (t '  )  dt '  /h. }  ;  ( 1-12 ) 
n  n  n  n 


the  resulting  coupled  equations  for  the  amplitudes,  an(t), 
are 


ifr  da . /dt  =  Z  v.P.  a  (t)  exp { — i ( e  -e . ) dt ' /h } ,  (1-13) 
y  n  jn  n  n  j 


. 


. 


where  we  have  used  the  relation 


-ifi  [9<f>  {r  ;R(t)  }/8t  ]->  =  -ifi  v.  (7)  +  <b  (r;R).  (1-14) 

.l  k  it  n 

Notice  that,  in  both  this  approach  and  that  of  equation 
(1-9),  transitions  arise  from  the  same  nonadiabatic  coup¬ 
ling  matrix,  P (R) .  The  connection  between  this  matrix  and 
the  adiabatic  criterion  can  be  seen  from  the  relation 

P  .  (R)  =  -iti  <<(>  .  I  (^  )  +  |  <p  > 

jn  T 3 1  R  r 1 yn 

=  -ih  «t>.  |  (^Dh  )-*|<j>  >/(e  -e.).  (1-15) 

rj'Rer|Tn'n  j 

Nonadiabatic  couplings  arise  from  the  changes  in,  or 
perturbations  of,  the  adiabatic  states,  as  the  nuclear 
configuration  alters  (hence  the  name,  "perturbed  station¬ 
ary  states"  theory).  The  presumption,  of  course,  is  that 
the  matrix,  P_(R)  ,  really  does  represent  the  effects  of 
physical  changes  in  the  wavef unctions .  However,  this 
presumption  is  not  correct,  as  I  shall  now  explain. 

3.  Electron  Translation  Factors 

Ih  order  to  understand  why  pss  theory  does  not 
provide  an  adequate  description  of  slow  collision  pro¬ 
cesses,  we  start  by  examining  the  asymptotic  (R-^00)  beha¬ 
viour  of  the  system.  Suppose  that  the  electron  is  in  a 
bound  atomic  state  on  nucleus  B,  asymptotically,  and  that 


* 


. 


II'  ^  ^ 
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this  "B-atom"  is  translating  with  respect  to  some  origin. 
Seen  from  this  origin,  the  electron  will  undergo  not  only 
its  internal  "atomic  orbital  motion,"  but  also  a  trans¬ 
lational  motion,  due  to  the  motion  of  the  nucleus  to  which 
it  is  bound.  Hence,  if  the  CMN  is  the  origin  of  electron 
co-ordinates  (as  in  the  molecular  system,  R,  r) ,  then  the 
wavefunction  must  include  a  factor  to  describe  the  trans¬ 
lational  motion  of  the  electron  with  respect  to  the  CMN. 

However,  pss  theory  gives  no  account  of  this  factor 
at  all:  the  Born-Oppenheimer  electronic  states  contain  no 
information  about  the  collision  velocity,  and  the  Xn(R) 's 
in  equation  (I-8a)  are  "nuclear  wavef unctions "  and  inde¬ 
pendent  of  the  electron  co-ordinates.  Effects  of  electron 
translation  must  somehow  be  buried  in  the  matrix,  P(R) , 
which  we  may  therefore  expect  to  show  some  unexpected  (and 
undesirable)  features. 

If  the  electron  is  in  an  atomic  state,  n,  on  nucleus 

B,  its  internal  motion  can  be  described  by  the  atomic 

orbital  ct  (r^)  ,  where  r„  is  the  vector  from  B  to  e  . 
rn  B  B 

Asymptotically,  this  orbital  is  totally  unperturbed  by  the 
distant  nucleus,  A,  and,  if  P(R)  is  to  represent  real  non- 
adiabatic  couplings,  we  would  therefore  expect  all  the 
matrix  elements  of  P(R)  to  tend  to  zero,  as  R+°°. 


In  fact,  P(R)  does  not  become  zero,  as  R-*°°;  moreover, 

the  non-zero  elements,  P  (°°)  ,  can  be  shown  to  arise  from 

B  -*■ 

simply  displacing  the  orbital,  4>n  (rB)  ,  with  respect  to 


■ 


. 


the  CMN ,  without  deforming  it  at  all.  These  displacement 
terms  are  wholly  spurious.  Unfortunately,  they  are  also 
often  large;  indeed,  in  some  cases,  they  dominate  even  at 
finite  R-values .  A  good  illustration  of  this  is  provided 
by  the  matrix  elements  of  P(R)  connecting  bound  molecular 
states  to  the  Born-Oppenheimer  electronic  continuum  (it 
can  be  shown^  that  such  couplings  are  related  to  the 
direct-impact  ionisation  cross  sections).  For  these 
transitions,  all  elements,  P^n(R) ,  go  to  zero  asymptotic¬ 
ally,  but,  in  the  interaction  region,  they  are  physically 

2 

unrealistic:  typically,  each  bound  state  is  coupled  to  30 

or  40  continuum  partial  waves,  and  the  envelope  of  coup¬ 
lings  extends  to  at  least  40  au.  One  of  the  aims  of  this 
thesis  is  to  show  that  a  careful  choice  of  the  "electron 
translation  factor"  eliminates  most  of  these  couplings 


entirely . 


A  further  problem  with  pss  theory  appears  if  we 
examine  P(R)  for  a  system  in  which  M^  /  Mg.  Because  the 
reference  origin  is  the  CMN,  P(R)  contains  terms  at  all  R- 
values,  which  depend  directly  on  the  nuclear  mass  asym¬ 
metry,  to  zero ' th  order  in  (m/p).  But,  in  actual  fact, 
the  electron  "sees"  only  the  electrostatic  potential  of 
the  nuclei  in  h  (r;R),  and  the  only  information  about  the 
nuclear  masses  that  can  appear  is  in  the  electronic 
reduced  mass,  as  a  correction  of  first  order  in  (m/p). 


This  anomaly  also  occurs  because  of  the  neglect  of 


'  riB  :  J  I  I  l  l  I  i 
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electron  translation,  and  has  been  discussed  in  great 
detail  for  HD+  by  Davis  and  Thorson.^^ 

Clearly,  the  matrix,  P_(R)  ,  does  not  represent  simply 
the  effects  of  real  change  in  the  electronic  wavef unctions , 
and  the  theory  needs  to  be  corrected.  To  see  how  this  can 
be  done,  we  examine  the  form  of  the  "electron  translation 
factor"  which  appears  in  the  wavef unction  asymptotically. 

If  the  electron  is  attached  to  B,  and  B  is  moving  at  a 
constant  velocity,  w_,  ,  with  respect  to  the  CMN,  then  (in 
the  time-dependent  language  of  equation  (1-5) )  the 
electron  translation  factor  (etf)  is  just 

Fb  =  exp{  (im/fi)  (wfi.r  -  h  wg  t)  } ,  (I-15a) 

and  the  asymptotic  electronic  wavefunction  is  not  simply 
<t>nB(rB),  but 

C  =  FB  *n  (?B>  •  (I-15b) 

The  factor,  FD ,  represents  just  the  momentum  and  kinetic 
energy  of  an  electron  moving  with  velocity,  w  ,  with 
respect  to  the  CMN.  Similarly,  if  the  electron  is  bound 
to  A,  there  is  a  corresponding  factor,  F  ,  and  velocity, 

"A. 

However,  in  the  interaction  region,  the  form  of  the 
etf  is  not  so  clear.  In  a  slow  collision,  we  want  to 


* 


describe  the  electron  motion  in  terms  of  the  Born- 
Oppenheimer  wavef unctions .  Now,  at  finite  R-values,  these 
states  are  "molecular,"  rather  than  "atomic,"  in  character, 
and  it  is  not  possible  to  say  that  the  electron  is 
attached  to  a  particular  nucleus.  Rather,  when  the  elec¬ 
tron  is  near  nucleus  A,  it  is  carried  along  by  A;  and, 
when  it  is  near  B,  it  translates  with  B.  Thus,  an  elec¬ 
tron  in  a  molecular  state  requires  a  "  molecular"  etf,  and 

g 

the  idea  of  a  switching  function  was  introduced  to  give 

the  electron  a  local  translational  velocity,  as  a  function 

7 

of  its  position.  Using  this  idea,  it  can  be  shown  that  a 
formally  correct  theory  of  slow  collisions  can  be  con¬ 
structed,  which  removes  all  the  defects  of  pss  theory, 
while  retaining  much  of  its  simplicity. 

4.  Switching  Functions  and  Translation  Factors 

The  rigorous,  fully  quantum  mechanical  derivation  of 

the  correct  coupled  equations  to  replace  (1-9)  is  very 
7b 

subtle.  The  defects  of  pss  theory  are  shown  to  arise 

•  •  •  ->■  ,  t 
from  an  inadequate  choice  of  co-ordinates,  (r,R) ,  in  which 

to  separate  the  internal,  "electronic"  motion  and  the 

translational,  "heavy-particle"  motion.  If  one  carefully 

selects  a  "scattering  co-ordinate,"  to  replace  R,  the 

(non-linear)  transformation  from  (r,R)  to  (r,£)  introduces 

important  new  terms  into  the  coupled  equations ,  terms 

which  remove  all  the  defects  of  pss  theory. 

In  this  thesis,  however,  I  prefer  to  use  an 


' 

. 
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alternative  derivation,  based  explicitly  on  the  etf 
7a 

concept,  which  brings  out  more  clearly  the  main  physical 

ideas.  I  will  show  that  the  etf  arises  naturally  from  the 

boundary  conditions  satisfied  by  Y  (or  T)  at  large  R,  and 

that  all  the  problems  with  pss  theory  can  be  solved  by 

attaching  an  appropriate  etf  to  each  Born-Oppenheimer 

state.  The  simplest  derivation  of  new  coupled  equations 

8  9 

assumes  the  classical  trajectory  approximation.  '  It 
will  turn  out  that  a  suitable  expansion  basis  is  given  by 
the  set 

|n(R)>  =  <f>n(r;R)  exp[(im/fi)  {v.  s-/t  ( 1-A2)  v2dt '  /  8 }  ]  , 

(1-16) 

where 


s  (r ; R)  =  H  {f(r;R)+A}  r 

y 


r  =  r  -  %  X  R 
9 


and  X  is  the  mass  asymmetry  parameter, 


(I-17a) 


(I-17b) 


A  =  ( ^A”Mb ) / ( Ma+Mq ) • 


(I-17c) 


f ( r ; R)  is  called  the  molecular  switching  function,  and  it 


allows  the  translational  motion  of  the  electron  to  vary 
with  its  position  in  the  molecule.  To  meet  the  boundary 


■ 
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conditions,  as  R+00 ,  f  must  approach  +1,  if  the  electron  is 
in  a  bound  state  of  B;  and  -1,  if  it  is  bound  to  A;  but  it 
is  otherwise  unconstrained.  The  question  of  how  f  should 
be  chosen  in  the  interaction  region  will  be  discussed  at 
length,  and  one  scheme  for  doing  so  is  the  main  topic  of 
this  thesis. 

Inclusion  of  the  etf  leads  to  important  modifications 
in  the  coupled  equations,  (1-9)  or  (1-13),  the  main  one 
being  the  replacement  of  the  matrix,  P(R) ,  by  the  cor¬ 
rected  nonadiabatic  coupling  matrix,  ?(R)+A(R),  where 

A.  (R)  =  ( im/fi )  <<J> .  |  [h  ,s  (r;R)  ]  |  <J>  >.  (1-18) 

J  *  A  I  C  X  A 


—V 

The  A-dependent  term  in  A(R)  exactly  cancels  the  spurious 
effects  of  mass  asymmetry  in  P(R) ,  and  transfers  the  origin 
of  electron  co-ordinates  to  the  geometric  centre  (r  is 
the  corresponding  co-ordinate).  In  the  limit,  R+°° ,  the 
term  containing  f  removes  the  spurious  asymptotic  coup¬ 
lings,  and,  in  effect,  refers  the  origin  of  co-ordinates 
to  the  nucleus  to  which  the  electron  is  actually  bound. 

5.  Ionising  Transitions 

In  the  earlier  discussion  of  the  adiabatic  criterion, 
I  pointed  out  that  electronic  continuum  states  cannot 
evolve  adiabat ically .  This  means  that  ionisation  cannot 
be  treated  conveniently  within  the  usual  framework  of  pss 
theory,  since  one  assumes  implicitly  that  the  basis 


■ 
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functions  behave  nearly  adiabatical ly .  However,  Thorson 
and  Levyi  have  shown  that  a  theory  of  ionisation  from 
adiabatic  bound  states  to  the  electronic  continuum  can  be 
developed  as  an  extension  of  the  pss  method.  Within  a 
first-order  perturbation,  or  distorted-wave,  approximation, 
they  proved  that  the  transition  probabilities  for  ion¬ 
isation  can  be  calculated  using  the  nonadiabatic  coupling 
matrix  elements  between  the  bound  electronic  states,  <f>  , 
and  the  Born-Oppenheimer  (fixed-nuclei)  continuum  states, 

(p  ( r;R ).  I  shall  give  a  very  brief  account  of  the  theory 
of  Thorson  and  Levy^  in  section  E  of  chapter  II. 

If  electron  translation  factors  are  ignored,  the 
matrix  elements  for  ionisation  are  just  the  couplings. 


P 

en 


(R)  , 


P 

en 


(R) 


(1-19) 


When  etf's  are  included  in  the  bound  state  description, 
these  couplings  are  replaced  by  P_en+Aen,  just  as  for 
couplings  between  bound  states. 

The  etf  corrections  are  of  critical  importance,  not 
only  to  the  calculation  of  ionisation  matrix  elements,  but 
also  in  more  general  slow  collision  problems;  in  partic¬ 
ular,  one  needs  a  way  of  choosing  the  switching  function 
in  the  interaction  region.  With  this  in  mind,  I  want  now 
to  describe  briefly  the  history  of  the  problem. 
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C.  History 

1.  Introduction 

In  common  with  so  many  other  problems  in  molecular 

phys  ics,  the  origins  of  slow  collision  theory  may  be 

14 

traced  to  the  famous  paper  by  Born  and  Oppenheimer  on 
the  quantum  mechanics  of  molecules.  Since  the  electron  to 
nuclear  mass  ratio  is  much  less  than  unity.  Born  and 
Oppenheimer  pointed  out  that  characteristic  nuclear  speeds 
in  the  bound  states  of  a  molecular  system  will  be  much 
lower  than  those  of  electrons;  therefore,  an  adiabatic 
separation  of  electronic  and  nuclear  motion  should  be  a 
good  approximation.  Their  approach  was  based  on  a  pertur- 

p 

bation  expansion  in  powers  of  (m/y)  4,  and  they  were  very 

careful  to  limit  their  discussion  to  non-degenerate, 

molecular  states,  in  which  the  nuclear  motion  is  bound; 

hence,  the  problem  of  asymptotic  boundary  conditions  never 

arose.  Nevertheless,  it  is  very  tempting  to  extend  their 

idea  directly  to  scattering  problems,  especially  since  the 

set  of  fixed-nuclei,  electronic  states  is  complete. ^ 

A  semiclassical  treatment  of  slow  collision  problems, 

12 

using  the  adiabatic  basis  states,  was  given  by  Mott 
(essentially,  equations  (1-12)  and  (1-13));  the  correspon¬ 
ding  quantum  mechanical  version  was  developed  by  Mott  and 
Massey11  (equations  (1-8)  and  (1-9))  and  named  the  per¬ 
turbed  stationary  states  theory.  In  1953,  Bates,  Massey 
1  f) 

and  Stewart  published  a  detailed  discussion  of  pss 
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theory,  applied  to  one-electron  problems.  As  far  as  I 

know,  this  paper  contains  the  first  recognition  that 

attempts  to  use  pss  theory  for  solving  slow  collision 

problems  will  run  into  fundamental  difficulties.  However, 

it  wasn't  until  1958  that  Bates  and  McCarroll  showed  that 

a  correct  theory  must  include  some  account  of  electron 

translation,  at  least  asymptotically. 

Since  then,  many  different  approaches  to  the  problem 

17 

have  been  tried.  In  1976,  Riera  and  Salin  wrote  an 
excellent  review  of  the  subject,  in  which  they  analysed 
the  various  approaches  and  their  defects,  and  pointed  out 
the  need  for  further  work.  Most  of  the  proposed  solutions 
and  all  of  them  which  correctly  meet  the  boundary  con¬ 
ditions,  may  be  seen,  in  retrospect,  as  corresponding  to 
different  choices  for  the  switching  function.  In  the  next 
subsection,  I  will  discuss  the  various  proposed  solutions 
from  this  viewpoint.  Of  course,  this  is  not  the  way  that 
people  actually  thought  about  the  problem,  but  it  is 
convenient  for  linking  the  development  of  the  subject  to 
the  present  work.  The  final  subsection  covers  the  recent 

7 

work  by  Thorson  and  Delos,  which  goes  a  long  way  towards 
presenting  a  unified  view  of  the  theory  of  slow  collisions 
2.  Proposed  Solutions  to  the  ETF  Problem 
Many  techniques  have  been  proposed  to  avoid  the 
difficulties  encountered  in  the  application  of  pss  theory 
to  practical  problems.  With  one  exception,  these  reduce 
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to:  (1)  techniques  which  take  f  =  constant;  these  are 
equivalent  to  using  a  single,  fixed  origin  for  electron 
co-ordinates,  and  (2)  methods  which  correspond  to  partic¬ 
ular  switching  functions.  In  the  first  category,  f  is  not 
a  switching  function,  as  it  cannot  meet  the  asymptotic 
constraints,  f-M-1  (-1)  on  B  (A)  ,  for  all  states;  only  the 
methods  in  the  second  category  are  capable  of  meeting  all 
boundary  conditions. 

a.  Methods  Using  a  Fixed  Reference  Origin 

i.  f  =  -  X.  For  completeness,  we  start  with 

unmodified  pss  theory,  which  is  recovered  if  we  take  f  = 

2 

-  X  (and  drop  the  term  in  ( 1—  X  )  in  equation  (1-16),  which 

7a 

is  just  a  phase) .  Thorson  and  Delos  point  out  that  pss 
theory  may  give  reasonably  good  answers  for  certain  types 
of  problems.  In  degeneracy  mediated  processes,  for 
example,  P^n(R)  may  be  very  large  in  the  region  of  near¬ 
degeneracy  (cf.  equation  (1-15)),  while  A^n(R)  is  much 
smaller,  and  it  may  be  an  adequate  approximation  to  ignore 
A ( R)  altogether. 

ii.  f  =  0.  This  has  the  effect  of  making  the  geo¬ 
metric  centre  (halfway  between  A  and  B)  the  origin  of 

9  18 

electron  co-ordinates,  and  Riley,  and  Schmid  have 
pointed  out  that  such  a  correction  may  often  be  useful. 

It  removes  the  spurious  effects  of  nuclear  mass  asymmetry 
from  pss  theory,  but  does  not  remove  the  remaining  errors, 
due  to  electron  translation  with  respect  to  the  geometric 
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iii.  f  =  +1.  This  choice  makes  nucleus  B  the  refer¬ 
ence  origin  for  electron  co-ordinates.  The  asymptotic 
wavef unctions  describing  electronic  states  bound  to  B  are 
then  exact,  but  those  associated  with  A  remain  incorrect. 

Matrix  elements  linking  two  states  on  B  have  the  correct 

17 

asymptotic  behaviour,  but,  as  Riera  and  Salin  point  out, 
there  is  no  improvement  for  couplings  involving  states  on 
A.  This  method  can  be  used  to  calculate  total  charge 
exchange  probabilities,  but  not  individual  transition 

cross  sections.  It  has  been  used  by  Piacentini  and 

19  20 

Salin,  and,  more  recently,  by  Winter  and  Lane,  to 

++ 

treat  charge  exchange  collisions  in  HeH 

21 

iVc  Scheme  of  Melius  and  Goddard.  These  authors 
attempt  to  correct  for  etf  effects  by  taking  an  origin 
which  is  effectively  a  weighted  average  over  a  molecular 
electronic  state  distribution.  While  ?(R)  will  generally 
vanish  asymptotically  for  electrically  asymmetric  systems, 
this  method  cannot  yield  a  consistently  correct  theory, 

and  it  will  certainly  fail  for  symmetric  systems,  such  as 

+  17  .  . 

H 2  .  Riera  and  Salin  have  given  a  detailed  criticism  of 

this  scheme. 

b.  Asymptotic  A  Approximation 

22 

Chen  et  al .  propose  to  correct  the  nonadiabatic 
couplings,  by  subtracting  from  P(R)  its  asymptotic  value, 
i.e.,  A  ( R)  =  A  (°°)  =  -  P(°°).  They  state  that  this  can  be 
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justified  by  a  formalism  similar  to  that  of  "renormali¬ 
sation"  in  quantum  electrodynamics;  it  does  not  correspond 
to  any  simple  assumption  about  f.  In  some  special  cases, 
this  may  be  a  useful  approximation,  but  it  is  not  general- 

ly  adequate,  as  is  pointed  out  both  by  Thorson  and  Delos , /a 

17 

and  by  Riera  and  Salin.  Many  pss  matrix  elements  are 
spuriously  large  at  finite  R-values,  even  though  they 
vanish  asymptotically,  and,  in  such  cases,  this  method 
offers  no  improvement  over  pss  theory.  The  ionising  trans¬ 
itions  provide  a  good  example  of  the  failure  of  this 
method. 

The  rest  of  this  subsection  deals  with  methods  in 
category  (2)  described  above,  those  which  satisfy  all 
asymptotic  boundary  conditions, 
c.  Asymptotic  f  Approximation 

This  is  the  method  introduced  by  Bates  and  Mc- 

3  . 

Carroll.  Equivalent  schemes,  expressed  in  a  variety  of 

formalisms,  have  been  developed  by  several  other  writers, 

2  3  24 

notably  Hahn  and  Matveyenko.  This  method  assigns  f  = 

+1  for  any  state  which  is  asymptotically  associated  with 

nucleus  B,  and  f  =  -1  for  states  associated  with  A.  For 

asymmetric  systems,  a  unique,  one-to-one  correlation 

exists  between  molecular  states  and  asymptotic  atomic 

states;  in  symmetric  systems,  however,  the  correlation  is 

3 

not  one-to-one.  Bates  and  McCarroll  proposed  a  modific¬ 
ation  to  handle  this  case,  but  it  leads  to  complicated 


■ 


'  * 


» 


non-orthogonalities  in  the  basis  states,  and  to  spurious 
couplings  between  states  of  £  and  u  symmetry. 

Even  for  asymmetric  systems,  at  finite  R,  this 

approach  does  not  give  a  very  satisfactory  description  of 

states  which  are  genuinely  "molecular"  in  character.  I 

will  compare  the  results  obtained  using  this  approximation 

with  those  from  the  more  sophisticated  switching  functions 

.  ++ 

reported  in  this  thesis.  For  HeH  ,  the  asymptotic  f 
approximation  is  very  good  in  some  cases,  but  poor  in 
others;  it  has  been  widely  used  in  actual  calculations, 
d.  Molecular  Switching  Function--Single  f 

In  order  to  accomodate  the  molecular  character  of  an 
electron  in  a  slow  collision,  Schneiderman  and  Russek^ 
introduced  the  idea  of  a  switching  function.  If  the  elec¬ 
tron  is  near  A,  f  =  -1;  and,  if  it  is  near  B,  f  =  +1.  To 
maintain  orthogonality,  they  used  a  single  switching 
function  for  all  electronic  states.  This  device  takes 
care  of  the  asymptotic  behaviour,  and  also  produces  a 
realistic  correction  in  the  interaction  region.  The  only 
problem  is  how  to  choose  f,  and  different  authors  have 
used  different  forms,  with  varying  success. 

g 

i.  Schneiderman  and  Russek 

2 

f  =  cos0/{l  +  (ot/ R)  },  ot  variable 


(1-20) 


I 


k '  . 


ii.  Mittleman  and  Tai 
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f  -  (1-S^)  { 1-exp ( -2an )}/{( 1+S2 )  [1+exp (-2an ) ] 

-  4  S  exp(-an)},  (I-21a) 


where 

S  =  {1  +  aR  +  (aR)2/3}  exp(-aR) ,  (I-21b) 

a  =  21/11,  D  =  (rA-rB)/R.  (I-21c) 

2  6 

iii.  Taulbjerg  et  al . 

f  =  (n/2)  (3-r]2)/{l  +  (a/R)2},  a  variable  (1-22) 

iv.  Levy  and  Thorson 

f  =  (rA2-rB2)/(rA2+rB2>  'X-23> 

f .  Molecular  Switching  Functions--state-dependent  f's 

If  we  are  prepared  to  push  the  physical  interpret¬ 
ation  of  the  switching  function  a  little  further,  then  we 
may  expect  f  to  be  dependent  on  the  bound  state  with  which 
it  is  associated.  As  I  pointed  out  earlier,  for  ionising 
transitions,  the  pss  matrix  elements  are  very  large  and 
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very  numerous;  the  corrected  couplings  are  also  very 

sensitive  to  the  choice  of  f,  and  this  was  one  of  the 

reasons  for  the  studies  of  switching  functions  made  by 

1  2 

Thorson  and  coworkers.  ' 

“t*  lb 

For  ,  Thorson  and  Levy  introduced  the  form 

f  =  tanh(B(R)  R  n),  (1-24) 

where  $(R)  is  to  be  chosen  for  each  initial  bound  state. 

Of  course,  if  we  use  state-dependent  f's,  the  electronic 

basis  is  once  again  non-orthogonal .  In  practice,  this 

2 

idea  was  remarkably  successful:  they  found  that  it  was 
possible  to  choose  3 (R)  such  that  the  number,  size,  and 
range  of  couplings  are  all  drastically  reduced.  Couplings 
to  all  but  the  first  two  or  three  continuum  partial  waves 
are  eliminated,  and  the  spurious  long-range  behaviour  of 
the  remaining  couplings  is  cancelled  (the  residual  coup¬ 
lings  have  a  range  of  about  10  au. ) .  The  corrected  matrix 
elements  to  the  higher  partial  waves  are  often  as  much  as 
four  orders  of  magnitude  smaller  than  the  corresponding 
pss  values. 

The  calculations  which  form  the  core  of  this  thesis 
confirm  these  results  and  extend  them  to  other  bound  states 
of  ,  and  to  systems  in  which  the  nuclear  charges  are 
unequal,  in  particular,  to  HeH++.  After  experimenting 
with  various  f's,  I  found  that  the  most  suitable  one  was 
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f  =  tanh{R[  (6/2)  {  (ZA+ZB)  f)+(ZA-ZB)  }  +  a  ln(ZB/ZA)]}, 

(1-25) 

where  a  and  3  are  variable  parameters.  The  method  used  to 

select  the  parameters,  a  and  3,  is  more  systematic  than 

2 

that  used  by  Thorson  et  al . ,  but  the  same  spectacular 
reductions  are  observed.  The  parameters  which  achieve 
this  are  independent  of  the  continuum  state  properties  and 
of  the  type  of  coupling  involved  (radial  or  angular) ,  but 
do  appear  to  be  characteristic  of  the  initial  bound  state. 
This  behaviour  is  illustrated  exhaustively  in  the  thesis, 
and  establishes  beyond  any  reasonable  doubt  that  the 
switching  function  approach  to  the  problem  of  electron 
translation  factors  has  a  sound  physical  basis. 

Since  the  switching  functions  determined  here  are, 
in  some  sense,  characteristic  of  their  corresponding  bound 
states,  it  is  reasonable  to  propose  their  use,  not  only  in 
ionisation  calculations,  but  also  in  more  general,  close¬ 
coupling  problems.  I  shall  discuss  this  proposal  at  some 
length  in  the  final  chapter  of  the  thesis. 

3.  Formal  Considerations 

One  difficulty  with  the  whole  subject  of  translation 

factor  corrections  has  been  the  lack  of  a  rigorous  and 

fully  quantum  mechanical  formulation;  most  discussions 

have  used  semiclassical  or  impact  parameter  treatments 

7 

only.  Recently,  however,  Thorson  and  Delos  have  pre¬ 
sented  a  formulation  which  runs  parallel  to  pss  theory. 
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remedies  all  its  defects,  and  establishes  the  correct  form 

for  the  coupled  equations  to  replace  (1-9).  They  obtained 

the  same  results  by  two  independent  derivations,  both  of 

which  use  a  single,  molecular  switching  function.  The 

first  method  introduces  an  etf  explicitly,  and  uses  a 

quantum  mechanical  transcription  formula  to  eliminate  the 

heavy  particle  velocity.  The  second  derivation7  develops 

27 

an  idea  originally  proposed  by  Mittleman:  the  switching 

function  is  used  to  define  a  non-linear  co-ordinate  trans- 
formation  from  molecular  co-ordinates,  (r,R),  to  new  co¬ 
ordinates,  (r,^) ,  and  the  etf  corrections  are  generated  in 

the  transformed  Hamiltonian.  In  effect,  this  approach 

14 

returns  to  the  original  idea  of  Born  and  Oppenheimer; 

the  new  "scattering  co-ordinate,"  %,  is  defined  in  just 

such  a  way  that  the  internal  electronic  motion  is  properly 

separated  from  the  translational  motion  of  the  electron 

and  nuclei,  and  a  separate  etf  is  not  needed. 

Some  aspects  of  the  problem  remain  unresolved. 

Other  than  the  asymptotic  constraints,  the  formulation  of 

7 

Thorson  and  Delos  gives  no  criteria  for  specifying  the 
switching  function.  Formally,  this  presents  no  difficulty; 
provided  one  uses  a  complete  set  of  electronic  states,  the 
solution  to  a  scattering  problem  must  be  independent  of 
the  choice  of  f.  However,  in  practice,  the  choice  of  f 
may  have  a  significant  effect  on  the  behaviour  of  P(R)  + 

A ( R) ,  as  the  results  presented  in  this  thesis  show.  It  is 


' 


. 


t 


.....  ■  -  '  ■  •  1  ”  i"*’ 


i  V  r 


not  known  whether  or  not  it  is  possible  to  give  a  unique, 
or  "optimum,"  definition  of  either  f  or  of  the  matrix, 

A(R) ,  at  finite  R-values.  Since  a  real  calculation  must 
use  a  severely  truncated  basis,  an  effective  scheme  for 
choosing  f  is  of  considerable  practical,  as  well  as 
formal,  importance. 

Using  the  classical  trajectory  approximation,  Riley 
2  8 

and  Green  have  used  the  Euler-Lagrange  method  to  obtain 
equations  which  should,  in  principle,  determine  the  form 
of  electron  translation  factors,  and,  hence,  f.  However, 
these  equations  are  very  complicated,  and,  so  far,  they 
have  only  been  solved  for  cases  where  f  does  not  depend  on 
the  electron  position. 

In  this  thesis,  I  have  taken  an  empirical  approach 

to  the  determination  of  f:  we  have  found,  by  calculation, 

that  it  is  possible  to  choose  an  f  for  each  bound  state, 

which  makes  a  large  number  of  couplings  from  that  state  to 

the  continuum  very  small.  We  think  that  the  existence  of 

these  "optimum"  f's  is  an  interesting  result,  and, 

possibly,  also  a  significant  one.  Since  these  f's  are 

state-dependent,  the  results  are  not  strictly  compatible 

7 

with  the  formulation  given  by  Thorson  and  Delos.  How¬ 
ever,  the  derivation  can  be  extended  to  include  the  case 
of  different  f's  for  different  states,  at  least  in  the 
classical  trajectory  limit;  this  will  be  covered  in  the 
Appendix. 
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D.  Reader's  Guide 


In  the  next  chapter,  the  basic  theory  of  electron 
translation  factors  is  described  in  detail.  After  he  has 
read  and  digested  this  chapter,  the  reader  is  advised  to 
return  to  chapter  I  and  re-read  section  B;  it  should  make 
a  lot  more  sense  the  second  time  around.  Chapter  II  is 
followed  by  two  highly  technical  chapters,  dealing  with 
the  computation  of  wavefunctions  (III) ,  and  matrix 
elements  (IV).  On  a  first  reading  of  chapter  III,  the 
reader  should  probably  look  at  sections  A,  B-l,  and  C-l 
only,  which  will  provide  useful  background,  without  too 
much  detail. 

Chapter  IV  is  a  different  matter.  Sections  A,  B, 
and  F  summarise  the  first  three  chapters  and  set  the  scene 
for  the  results.  I  am  very  much  afraid  that  sections  C,  D 
and  E  will  defeat  all  but  the  most  dedicated  reader,  since 
they  deal  in  great  detail  with  the  computation  of  indiv¬ 
idual  matrix  elements,  and  hence  consist  mostly  of  algebra 
They  provide  a  recipe  for  anyone  wishing  to  carry  out 
similar  calculations,  but  should  be  read  (if  at  all)  only 
when  the  reader  feels  he  has  understood  the  rest  of  the 
thesis . 

The  results  of  my  work  are  presented  in  chapter  V. 
The  way  in  which  the  switching  functions  were  chosen  is 
described,  and  some  pictures  of  these  "optimum"  f's  are 
given.  These  are  followed  by  some  two-dozen  pictures  of 
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the  resulting  "optimised"  matrix  elements.  It  is  ultim¬ 
ately  these  pictures  which  provide  the  justification  for 
the  etf  approach  to  slow  collision  theory.  The  systematic 
cancellation  of  the  pss  matrix  elements  by  several  orders 
of  magnitude  is  not  fortuitous.  In  the  final  section  of 
this  chapter,  some  of  the  implications  of  these  results 


are  discussed. 
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CHAPTER  II 


THEORY 

A.  Co-ordinates 

Consider  a  system  composed  of  three  particles,  with 
masses,  nu  ,  i  =  1,  2,  3;  and  co-ordinates,  R^ ,  with 
respect  to  an  external  origin,  and  r^ ,  with  respect  to  the 
centre  of  mass  of  the  system.  Some  other  co-ordinates 
will  also  be  useful  (see  figure  1) . 
i.  Definition  of  the  centre  of  mass: 


E  .  m.  r .  =  0 . 

ill 


,  .  -»■  0 
li.  Centre  of  mass,  R^..: 

CM 


*CM  =  2i  mi  V/'V 


-*■  0  +  0 
ri  "  Ri  -  rcm' 


=  E  .  m .  . 

l  l 


■*  0 

iii.  Centre  of  mass  of  A  and  B,  R 


R^  =  (m1R1°  +  m2R2°  )  /  (m1+m2 )  . 


(II-l) 


(II-2) 

( I 1-3) 


( II-4) 
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Figure  1.  System  Co-ordinates 


(II-5a) 


iv.  Relative  co-ordinates,  q, ,  q  : 

1  2 

ql  =  r2  "  rl' 

q2  =  r3  (mlrl  +  m2r2)/(mi+m2) •  (II-5b) 

It  is  well  known  that  we  can  write  the  total  kineti 
energy  of  the  system  as 


T 


(II-6a) 


where 


TCM  =  35  It1T(RCM)2'  (II-6b) 

T'  =  %  y1(q1)2  +  %  y2(q2)2,  (II-6c) 

yl  =  m1m2/(m1+m2)  ,  y2  =  m3  (m^n^) /m^, .  (II-6d) 


We  can  therefore  factor  out  the  centre  of  mass  motion  and 
write 


L  = 


y(5l'52>  exp(iKc°.Rc°) 


( II-7) 


The  internal  motion  is  described  by  the  co-ordinates,  q^, 
the  vector  joining  two  of  the  particles;  and  q2 ,  the 
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position  vector  of  the  third  particle,  measured  from  the 
of  mass  of  the  other  two. 

Clearly,  there  are  three  ways  in  which  these  co¬ 
ordinates  can  be  chosen.  For  systems  composed  of  two 
heavy  particles,  masses  ma.  Mg,  and  an  electron,  mass  m  , 
they  are  called  the  "molecular,"  "A-atom,"  and  "B-atom" 
co-ordinates.  These  are  illustrated  in  figure  2. 
i.  Molecular  co-ordinates,  R,  r  (figure  2a): 

a  =  r  + 

ql  R'  =  r;  (II- 8a) 

y  =  MaV(MA+MB}'  m  =  m0  (MA+MB)/MT;  ( II-8b) 

MT=MA+MB+m0'  X=  (Ma“MB)/(MA+MB)*  (II“9) 

A  is  called  the  mass  asymmetry  parameter.  It  is  also  use¬ 
ful  to  define  r  ,  the  position  vector  of  the  electron, 
measured  from  the  geometric  centre  (halfway  between  A  and 
B)  s 

Tg  ~  r  ~  %  A  R.  (11-10) 

ii.  A-atom  co-ordinates,  R  ,  r  (figure  2b): 

=  =  ^MAMT/ <MA+tV  "  fm0/(MA+m0)  t11-113) 


nt  '  -  :l  °~ 
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(c) 


Figure  2.  Jacobi  Co-ordinates 


(II-llb) 


-  r  -  %  (X-l)  R; 

PA  =  MB(MA+mo)/MT»  mA  =  MAm0/(MA+m0)  .  (11-12) 

iii.  B-atom  co-ordinates,  ^ ,  ?B  (figure  2c): 

q2  =  RB  =  ^MBMT/(Ma+MB)  }R  +  {m0/(MB+m0)  >?»  (II-13a) 

ql  =  rB  =  r  ”  **  (*+1>  R;  (IX-13b) 

PB  =  MA(MB+rn0)/MT'  mB  =  MBm0/(MB+m0)  •  f11"14) 

I  want  now  to  introduce  a  function,  f  ,  where  J  is  A 
or  B,  defined  such  that  ffi  =  +1,  and  fft  =  -1.  with  this 
device,  we  can  use  a  single  form  to  express  both  sets  of 
atomic  co-ordinates,  in  terms  of  molecular  co-ordinates. 


RJ  =  <H/Kj> 

{R  +  (m/2y) (fj+X)  ?}, 

(II-15a) 

. 

tj  =  r  -  Js 

(fj+X)  R; 

(II-15b) 

Mj/v  =  {1  + 

(m0/2u)  (1+fjX)  )  (Ma+Mb)/mt, 

(II-16a) 

Vmj  -  1  + 

(mQ/2p)  (1+fjX)  . 

(II-16b) 

Equations  (IX  15)  and  (11-16)  airc  nssdsd  foir  the  naxt 
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section,  in  which  I  examine  the  boundary  conditions  on  ¥ 

and  show  why  the  pss  theory  gives  a  poor  description  of 

collision  processes.  Two  other  equations  will  also  be 
useful: 


V(W  =  1  +  *  m0/M, 


(II-17a) 


M 


aMb/Mt  =  =  y. 


(II-17b) 


B *  Asymptotic  Boundary  Conditions 

and  Translation  Factors' 

1.  Origin  of  the  Translation  Factor73 
Suppose  that  asymptotically  the  electron  is  in  the 
atomic  state,  <f> n  ,  on  atom  J,  with  energy,  e  .  Then  the 

J  n  j 

asymptotic  solution  to  the  time-independent  Schr&dinger 
equation  is  ' 


0 


¥  -v  <(>„  (rj)  expfik  .R  )  , 


(11-18) 


where 


2u  (El-e 
J  n 


(11-19) 


We  now  transform  the  plane  wave  factor  to  molecular  co¬ 
ordinates,  using  equation  (II-15a) .  However,  this  must  be 
done  carefully,  since  the  appearance  of  the  atomic  reduced 
j,  in  equation  (11-19)  means  that  the 


mass,  y 


wave-vector , 


.  .  ■  .  ••  ’• 


k j ,  is  an  intrinsic  function  of  the  atomic  configuration. 

(The  extrinsic  dependence  on  e  is  not  important.)  This 

J 

configuration  dependence  can  be  eliminated  by  defining  a 
molecular  wave-vector,  k,  such  that 


(11-20) 


and  this  implies  that 

kj*Rj  =  k*{R  +  (m/2y)  (fj+X)  r}.  (11-21) 

But  (m/y)  <<  1,  so  we  expand  (y/y  )  in  powers  of  (m/y)  , 
and  keep  only  the  lowest  terms.  From  equations  (II-16a) 
and  (II-17a) ,  it  can  be  seen  that 

(U/Hjp  =  1  -  (m0/8y)  (l+2fJA+A2),  (11-22) 


and 


kj.Rj  -  k.{R  +  (m/y)  s}, 


(11-23) 


where 


s  =  h  (f  +A)  r  -  ( 1 — A 2 )  R/8. 

u  y 


(11-24) 


The  term,  exp  [  (im/u)  £  . s  ]  ,  which  thus  appears  in 


■ 
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equation  (11-18),  is  called  the  electron  translation  fac¬ 
tor;  and,  if  we  define  a  classical  heavy-particle  velocity, 
■+  -> 

v  -  hk/p,  it  can  be  given  a  simple  physical  interpretation. 
(fJ+A)v/2  is  just  the  translational  velocity  of  the  elec¬ 
tron,  with  respect  to  the  centre  of  mass  of  the  nuclei, 
due  to  the  motion  of  the  nucleus  to  which  it  is  bound;  and 
the  remaining  term  represents  the  electron  transport 

kinetic  energy,  due  to  the  nuclear  motion.  Noting  that 
2 

fj  -  1,  for  a  straight-line  trajectory  R  =  b  +  vt  (b  is 
called  the  impact  parameter) ,  we  can  re-write  equation 
(11-23)  as 

kJ'Rj  ~  k,R  +  (m/k)  (w.r  -  w2t/2) ,  (II-25a) 


where 


w  (fj+AJv/2.  (II-25b ) 

Part  of  the  kinetic  energy  is  independent  of  the  system 
configuration  (i.e.,  of  fj) ,  and  it  is  convenient  to 
isolate  this  term: 


’V’Ar  =  [R+  (m/u)s’  ]  -  (m/ti)  (1-A2)  v2t/8,  (II-26a) 


where 


. 


. 
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S'  =  *5  (fT+X)  r  . 

J  g 


( II -2 6b) 


Alternatively,  we  could  absorb  the  final  term  into 
the  definition  of  the  molecular  wave-vector.  If  we  take 


*.T  =  k', 


(II-27a) 


where  jj  is  defined  in  equation  (II-17b)  ,  and  notice  that 


mQ/y  =  m/y, 

then  it  is  easy  to  see  that 

kj • Rj  -  k ' . [R  +  (m/y ) s']. 


(II-27b) 


(11-28) 


2.  Defects  of  PSS  Theory 
a.  Introduction 

I  want  now  to  examine  the  pss  theory  in  the  light  of 
this  study  of  the  boundary  conditions.  As  the  reader  will 
recall,  in  the  pss  method,  the  wavefunction  is  expanded  in 
terms  of  the  Born-Oppenheimer,  molecular  eigenfunctions 
(equation  (1-8)): 11 

f(r,R)  =  £nxn(R)  <Pn(r;R)  ,  (11-29) 

and  the  resulting  equations,  (1-9)  ,  are  coupled  by  the 
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matrix,  P(R),  where 


(R) 


(11-30) 


In  pss  theory,  individual  terms  in  the  expansion, 
(11-29),  are  interpreted  asymptotically  as  scattering 
wavefunctions  for  electronic  states,  <(>n"  [i.e.,  the  states 

*n 'r'R'R"*’0°n  »  and  these  states  are  presumed  to  correspond 
in  some  simple  way  to  the  exact  atomic  states,  4  For 

^  J 

simplicity,  let  us  first  consider  electrically  asymmetric 
systems,  in  which  this  correspondence  is  one-to-one.  If 
one  can  express  the  exact  plane-wave  state,  (11-18),  as  a 
single  product,  <f>n  Xn(R),  then,  and  only  then,  one  can 
interpret  this  product  as  a  scattering  wavefunction  for 
the  atomic  state,  d> 

nj 

However,  it  is  immediately  obvious  that  no  electron 

translation  factor  appears  in  the  product  <p  \  (R)  .  For 

our  purposes,  this  is  the  important  problem,  but  there  is 

another,  minor  error,  related  to  the  fact  that  d>  °°  is  not 

n 

exactly  equal  to  <J>n  ;  for  completeness,  we  will  dispose  of 
this  latter  defect  first, 
b •  Asymptotic  Electronic  States 

The  exact  atomic  states  are  eigenfunctions  of  a 
Hamiltonian,  h_^  ,  which  uses  the  atomic  electron  reduced 
mass,  m  ,  whereas  the  molecular  reduced  mass,  m,  appears 
in  the  Hamiltonian,  h  . 


Thus,  asymptotically. 


. 


. 
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hat  _  he  _  (ft2/2)  ("ij-1  -  ra  l)  (Vr2)g; 


(11-31) 


and  the  final  terra  can  be  treated  as  a  small  perturbation, 
as  it  is  of  order  (m/y)  times  the  electron  kinetic  energy. 
It  then  follows  that  the  atomic  eigenvalue,  en  ,  is  given, 
correct  to  first  order,  by  ^ 


£nj  “  {£n(R)  "  tt2/2)  (nrij'^m  X)  <4>n  I  (Vr2 )  +  |  <f>n > } , 

(II-32a) 

where  <(.  tends  to  the  approximate  atomic  state,  6  “ 

n  * 

Using  equations  (11-15)  and  (16),  we  can  rewrite  the 
second  term,  so  that 


lim 


n 


{en(R)  ‘  /2p)<4>n|  (VR  )jUn>}. 


(II-32b) 


As  equation  (1-11)  shows,  this  correction  does  appear  in 
the  pss  coupled  equations,  (1-9),  a  fact  which  is  well- 
known;  as  one  might  expect,  it  does  not  appear  in  the 

classical  trajectory  equations,  (1-13).  Unless  one  is 
interested  in  isotope  splitting  effects,  these  corrections 
to  the  atomic  binding  energies  are  unimportant;  I  propose 
to  ignore  them  and  assume  that  the  molecular  states  pro¬ 
vide  an  adequate  asymptotic  description  of  the  limiting 

atomic  states.  In  their  rigorous  formulation,  Thorson  and 
7b 

Delos  treat  this  point  much  more  carefully,  and  show 
that  the  corrected  theory  also  recovers  the  correct 
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binding  energies,  to  the  required  accuracy  in  (m/y) . 

In  electrically  symmetric  systems,  the  Born- 
Oppenheimer  wavef unctions  occur  in  asymptotically  degen¬ 
erate,  (g,  u)  pairs.  However,  limiting  atomic  states  can 
be  constructed  as  linear  combinations  of  the  g  and  u  mol¬ 
ecular  states, 

‘t’rij  =  2  ^  f<t>ng(?;R)  +  fj  *nU(?;R)},  (11-33) 

afber  which,  the  discussion  proceeds  as  before. 
c*  Translation  Factors  and  PSS  Theory 

Let  us  now  turn  our  attention  to  the  really  impor¬ 
tant  defect  of  pss  theory.  Suppose  that  the  electron  is 

asymptotically  in  state  d>  ,  where 

nj 

4>n(r;R)  -v.  4>n  (?j)  ,  as  R-«».  (11-34) 

J 

An  asymptotic  solution  to  the  pss  equations  (1-9)  may  be 
written 

y  ^  <J)n  exp(ik.R)  ,  (11-35) 

J 


where 


1a2k2 


2y  [E-en  (°°)  ]  ; 


(11-36) 
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..  . 
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and,  in  pss  theory,  this  is  assumed  to  be  equivalent  to 
(H-18).  However,  from  (11-26),  or  (11-28),  it  can  be 
seen  that  no  translation  factor  appears,  and  hence,  Y  as 
defined  by  (11-35)  is  not  equivalent  to  (11-18) .  Only  by 
using  a  complete  expansion  (including  the  electronic 
continuum) ,  can  one  make  the  form,  (11-29) ,  represent  a 
solution  like  (11-18) ;  such  an  expansion  is  clearly  im¬ 
practical.  Since  individual  terms  in  the  pss  expansion  do 
not  represent  real,  physical  states,  it  is  not  surprising 
that  the  couplings  between  them  show  unphysical  behaviour. 

To  show  how  the  spurious  asymptotic  coupling  arises, 
let  us  examine  the  nonadiabatic  coupling  matrix,  P(R),  in 
detail.  P(R)  is  the  matrix  of  the  operator,  -iti(\?  )  +  ,  in 
the  space  of  the  Born-Oppenheimer  functions,  and  it  is 
evaluated  holding  the  electron  fixed,  with  respect  to  the 

centre  of  mass  of  the  nuclei.  From  (II-15b) ,  it  can  be 
seen  that 

"lft(Vr  =  f(VR}?  “  *  (fj+A)  (Vr  )gl.  (11-37) 

J  j 

Asymptotically,  the  contribution  to  P(R)  arising  from  the 
first  term  in  equation  (11-37)  represents  deformation  of 
the  adiabatic  wavefunctions ,  as  R  changes,  and  it  must 
vanish,  since,  according  to  equation  (11-34),  the  <J>  (r;R) 
become  functions  only  of  r^..  However,  some  states  asymp¬ 
totically  associated  with  the  same  centre  will  be 


■ 


' 


..  ^ 


connected  by  the  momentum  operator,  )  +  ,  and,  for 

r  J  ^ 

such  states,  Pjn(R)  will  tend  to  a  constant,  as  r-h».  This 
term  represents  the  displacement  of  a  wave function ,  with¬ 
out  deformation,  and  it  is  entirely  spurious.  The 
physically-appropriate  reference  origin  is  atom  J,  not  the 
CMN.  At  finite  R-values ,  it  is  not  possible  to  decompose 
— (R)  uniquely  into  deformation  terms  plus  displacement 
effects,  but  it  is  clear  nevertheless,  that  the  latter  are 
often  large,  and  that  they  have  no  place  in  the  theory. 

3.  Modified  Electronic  Basis 

From  the  above  discussion,  it  is  clear  that  we  can 
avoid  the  defects  of  pss  theory  by  including  an  account  of 
electron  translation  in  the  electronic  basis.  In  the 
classical  trajectory  approximation,  we  can  do  this  simply 

by  replacing  the  Born-Oppenheimer  basis  with  the  modified 
states , 


■tn{?;R(t)}  =  0n(?;R(t)}  Fn,  (11-38) 

where  Fn  is  an  electron  translation  factor,  defined  by 

Fn  =  exp {(im/h)  v.sn>.  (11-39) 

The  exponent,  s^,  must  tend  to  the  correct  asymptotic 
value,  given  by  equation  (11-24) ,  but  it  should  also  re¬ 
flect  the  fact  that  the  electron  is  molecular  for  much  of 
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the  collision,  with  appreciable  density  near  both  nuclei. 
For  this  reason,  we  introduce  a  molecular  switching 
function,  f (r;R) , ^  defined  such  that 


lim 

r^oo  f(r;R)  =  fj,  finite; 
r^oo  f(r;R)  =  0,  R  finite; 


(II-40a) 


(II-40b) 


f(-r;R)  =  -  f (r ; R) ,  if  Z  =  Z 


B 


(II-40c) 


but  f  is  otherwise  unconstrained.  We  will  thus  obtain  a 
suitable  etf  by  defining  s  ,  such  that 


V*sn  =  v’s  "  ^  (H  )v^dt'/8, 


(II-4 la) 


s  =  H  (f+X)  ?  . 

g 


(II-4 lb) 


Note  that  the  s  of  equations  (11-41)  tends,  asymptotically, 

to  the  s'  of  equation  (11-28).  In  equation  (11-38),  I 

have  allowed  for  different  Fn ' s  for  different  states;  this 

flexibility  would  be  provided  by  choosing  distinct  f  ' s  in 

n 

the  definition  of  s.  However,  for  the  moment,  I  shall 
take  the  simpler  case  of  a  single,  common  f,  and  genera¬ 
lise  the  results  obtained  later. 

From  the  point  of  view  of  a  fully  quantum  mechanical 
theory,  the  modified  basis  is  still  defective,  since  the 
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classical  nuclear  velocity  appears  in  the  translation  fac¬ 
tor.  In  their  first  formulation,  Thorson  and  Delos 7a 
introduce  an  ad  hoc,  quantum  mechanical  transcription  to 
eliminate  v;  their  second  derivation,713  based  on  a  co¬ 
ordinate  transformation,  gives  the  same  final  result  and 
justifies  the  transcription.  Here,  I  shall  use  a  classical 
trajectory  approach,  which  is  much  simpler,  but  still 
retains  all  the  essential  features. 

C.  Classical  Trajectory  Equations 
this  section,  I  will  derive  the  appropriate 
coupled  equations  to  replace  equations  (1-13) ,  using  a 
basis  of  Born-Oppenheimer  adiabatic  states,  modified  by  a 
translation  factor.  We  assume  that  the  nuclei  follow  a 
specified  classical  path,  R(t) ,  with  velocity  v  =  dR/dt, 

and  write  the  solution  to  the  time-dependent  Schrddinger 
equation  as 

T  =  Zn  an(t)  exp{-i/ten  (t ' )  dt '  /h } ,  (n-42) 

where  $n(r;R)  is  defined  by  equations  (11-38)  to  (11-41). 
This  implies  the  system  of  coupled  equations 

ifi  b  =  [h  +  v.D]  b_,  (11-43) 


where 


■ 


m 


(11-44) 


bn  =  an  exp{-i/ten(t')dt,/fi}, 


h  .  =  <4>  .  |h  I  $  > , 

jn  j 1  e 1  n 


(11-45) 


V  = 


(11-46) 


Now 


h .  —  e  6.  +  <<f>  .  I  f  [h  » F 1  ld>  > 

jn  n  jn  1  1  g  ' r  *  l  * 


(11-47) 


P4  n  (R)  =  -ifc  <(p  .  I  (  V)  ■>  I  <f>  >; 
jn  j  1  R  r 1 


(11-48) 


and  the  remaining  term  is  thus 


<4>j|F  {[he,F]  -  ifi  9F/  3t }  [  4>n>  . 


(11-49) 


If  we  discard  the  terms  in  (11-49)  involving  9v/9t,  and 

2 

those  of  order  v  ,  we  find  that 


ifi  b  =  [e  +  v.  (p+A)  ]  b 


(II-50a) 


or 


ift  a.  -  En  v.(Pjn+Ajn)  an(t)  exp  (en-e  . )  dt  -  /k  } , 


(II-50b) 


where 
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Ajn(R)  -  (im/h)  <4>j  |  [he,s]  |  t^n>.  (11-51) 

Equations  (11-50)  are  the  desired  classical  trajectory 
equations.  The  only  difference  between  (11-50)  and  the 
pss  equations,  (1-13),  is  that  the  nonadiabatic  coupling 
matrix,  P(R) ,  has  been  replaced  by  P(R)+A(R) 

D.  Effect  of  Translation  Factor 
Corrections 

From  equation  (11-51) ,  it  can  be  seen  that 

Ajn(R)  =  (ej_en)  <(^jl^l^n>  (II-52a) 

=  -iti  <^j|(Vrs).Vr  +  k  (Vr2s)  |  0n>  ;  ( II-52b) 

and  the  reader  will  recall  that 


s  =  h  (f+A)  r 


(11-53) 


Therefore,  as  R-*-°°, 


A.  (R)  (f,+  X)  <<(> .  I  )U  > 

311  J  V  rj  nJ 


(II-54a) 


But,  from  equation  (11-37) ,  we  see  that 


p,n(R)  ^  -  H  (fT+X)  <<f>,  |  (-ihv  )U  >. 

J  J  J  L  J  nJ 


(II-54b) 
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Thus,  A (R)  exactly  cancels  the  spurious,  infinite-range 
couplings  contained  in  P(R): 

r+oo  +  A(R)  ]  =  0.  (11-55) 

Of  course,  that  was  the  main  reason  for  introducing  the 
etf  in  the  first  place! 

For  systems  with  unequal  nuclear  masses,  P(R)  con¬ 
tains  a  term  proportional  to  X,  since  the  origin  of  elec¬ 
tron  co-ordinates  is  the  centre  of  mass  of  the  nuclei. 
However,  nonadiabatic  couplings  are  supposed  to  represent 
properties  of  the  electronic  wavef unctions  and  thus  should 
be  independent  of  X.  The  additional  term  is  fictitious, 
and  it  is  exactly  cancelled  at  all  R  by  a  corresponding 
term  in  A(R) .  This  is  easily  shown: 

=  (~ifiVR)j  -  H  x  (-ifiV  )r,  (II-56a) 

g  x 


and 


(im/fr)  [h  ,  s  ]  —  (im/fi)  [h  ,s  ]  +  ^X(—  ifW  );>,  (II  — 56b) 

t;  g  r  K 


where 


s 


g 


h  f  r  . 
g 


( II-56c) 
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The  second  terra  in  (II-56b)  exactly  cancels  the  second 
term  in  (Il-56a) ,  so  we  have  the  result  that 


P(R)  +  A(R)  =  Pg(R)  +  Ag(R), 


(II-57a) 


where 


(R) 


=  -ifi  <<f>  .  |  ($)-» 

J  1  R  r 

g 


(II-57b) 


and 


Ajn(R)  "  <0j  Uhe,sg]  |<J>n>.  (II-5  7c) 

How  is  A (R)  to  be  interpreted?  At  large  R,  we  have 
seen  that  P(R)  can  be  separated  into  a  distortion  term, 
which  tends  to  zero,  and  a  displacement  tern.  The  latter 
is  spurious,  and  it  is  exactly  cancelled  by  A(R) .  We  may 
suppose  that  A(R)  plays  much  the  same  role  at  finite  R- 
values:  it  identifies  and  removes  that  part  of  P(R)  which 
represents  displacement  without  deformation,  leaving  only 
the  distortion  effects.  However,  since  there  is  no 
obvious  and  unique  way  of  isolating  just  the  translation 
effects ,  A ( R)  can  only  be  defined  uniquely  in  the  asymp¬ 
totic  limit.  This  ambiguity  is  reflected  in  the  fact  that 
the  theory  gives  no  way  of  determining  f  in  the  inter¬ 
action  region.  Formally,  this  does  not  matter:  the  set  of 
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Born-Oppenheimer  states  is  complete,  and  any  "sensible"  f 
which  satisfies  the  boundary  conditions  will  do. 

In  practice,  the  choice  of  f  is  important,  since  one 
must  always  work  with  a  truncated  basis,  and,  in  many 
cases,  the  corrected  matrix  elements  are  very  sensitive  to 
the  choice  of  f.  In  this  work,  I  examine  the  way  in  which 
the  corrected  matrix  elements  for  ionising  transitions 
depend  on  the  switching  function.  The  behaviour  of  these 
matrix  elements  allows  us  to  choose  very  clear-cut, 
"sensible"  f's. 

Two  tasks  are  required  to  complete  this  chapter.  In 
section  E,  I  give  a  brief  account  of  the  theory  of  direct- 
impact  ionisation  in  slow  collisions,  to  provide  a  back¬ 
ground  for  the  examination  of  ionisation  matrix  elements. 
Finally,  section  F  deals  with  the  transformation  of  elec¬ 
tron  co-ordinates  from  a  space-fixed  reference  frame,  to  a 
frame  rotating  with  the  molecular  axis. 

E.  Ionisation  in  near-adiabatic 

Collisions 

1.  Introduction 

As  I  pointed  out  in  chapter  I ,  ionising  transitions 
cannot  be  handled  within  the  close— coupled  expansion 
scheme  of  pss  theory  (even  if  electron  translation  effects 
are  included) .  The  expansion  uses  a  truncated  basis  of 
adiabatic  states,  and  assumes  that  (in  a  slow  collision) 
electrons  in  those  states  behave  nearly  adiabatically . 
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But  r  because  of  their  high  degeneracy,  the  continuum  and 
Rydberg  levels  can  never  behave  adiabatically ,  and  to 
include  them  in  a  straightforward  pss  expansion  would  mean 
using  an  intractably  large  basis. 

However,  in  direct  impact  processes,  large  changes 
in  the  electronic  energy  occur  upon  transition,  and  hence 
the  probability  of  such  transitions  is  small,  at  least  in 
slow  collisions.  Thus,  one  may  be  able  to  treat  ionis¬ 
ation  from  tightly  bound  electronic  states  by  a  pertur¬ 
bative  approach,  in  which  the  perturbation  causing  trans¬ 
itions  is  the  "deformation"  of  the  initial-state  wave- 
function,  projected  onto  the  continuum  states  of  the 
system.  What  is  then  needed  is  a  description  of  these 
states,  and  a  definition  of  the  probabilities  of  trans¬ 
itions  to  them. 

A  first-order  perturbation  theory  of  direct  impact 
ionisation  in  slow  atomic  collisions  has  been  given  by 
Thorson  and  Levy . It  leads  to  a  "distorted-wave" 
approximation  for  the  ionisation  cross  sections,  and  the 
reader  is  referred  to  their  long  paper  for  a  detailed 
discussion;  here,  I  shall  describe  only  the  main  physical 
ideas . 

2.  The  "Fast  Electron"  Approximation 

An  adiabatic  electron  is  one  whose  (bound)  orbital 
speed  is  so  much  greater  than  the  nuclear  speed  that  it 
can  adjust  smoothly  to  changes  in  the  nuclear  position. 
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By  contrast,  a  continuum  electron  has  a  classically 
P  iodic  orbit,  so  that  it  cannot  be  said  to  follow  the 
nuclear  motion  adiabatically .  Nevertheless,  the  charac- 
ristic  speed  of  a  continuum  electron  will  still,  in 
general,  be  much  greater  than  that  of  the  nuclei.  This 
means  that  an  electron  in  an  unbound  orbit  will  move 
through  the  system  and  out  to  infinity  in  a  time  so  short, 
that  the  nuclear  configuration  will  hardly  change  at  all, 
and  the  orbit  will  essentially  be  that  of  an  electron  in 
the  fixed— nuclei "  system. 

In  the  language  of  quantum  mechanics,  an  electron 
stantaneously  excited  to  a  continuum  state  at  time,  t', 
will  produce  a  signal  at  a  detector  at  some  later  time,  t; 
but  the  Signal  will  reflect  the  nuclear  configuration  at 
the  excitation  time,  5(f),  and  not  that  at  the  signal 
time,  R(t).  By  taking  a  coherent  sum  of  the  instantaneous 
amplitudes  at  each  point  along  a  trajectory,  R(t'),  one 
obtains  the  cumulative  signal  amplitude  over  an  entire 
collision.  The  instantaneous  amplitudes  for  the  excita¬ 
tion  of  an  electron  to  a  continuum  state  with  energy,  e , 
is  the  projection  of  the  nonadiabatic  perturbation  of  the 
initial  state  onto  the  Born-Oppenheimer  continuum  state; 
and  the  matrix  element  for  this  is  just  P  (in  pss 
theory) ,  or  P£n  +  A£n  (in  the  corrected  theory) .  To  give 
a  correct  description  of  the  ionised  electron,  as  it  moves 
out  to  the  detector,  one  must  also  include  a  transformation 
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from  the  rotating  molecular  frame  back  to  the  laboratory 
frame,  but  the  basic  physical  idea  of  the  Thorson  and  Levy 
ionisation  theoryla  is  that  described  above. 

On  either  side  of  the  ionisation  threshold,  there  is 
a  band  of  states  in  which  the  electron  is  neither  "fast" 
nor  adiabatic."  In  these  states,  characteristic  electron 
speeds  are  of  the  same  order  as,  or  less  than,  the  nuclear 
speeds,  and  hence,  for  collisions  involving  these  states, 
the  electron  motion  will  be  very  complicated.  Second 
order  processes,  such  as  ionisation  followed  by  recapture, 
will  be  common,  and  no  simple  theory  exists  to  deal  with 
these  "slow  electron"  levels.  For  slow  collisions,  the 
band  of  states  is  less  than  0 . 5  eV  wide,  and  transitions 
to  such  states  are  unlikely  to  contribute  much  to  the 
overall  cross  section. ^a 

3.  Ionising  Transitions  and  Switching  Functions 

For  the  purposes  of  this  thesis,  the  important  con¬ 
clusion  to  be  drawn  from  the  work  of  Thorson  and  Levy^a  is 

that  the  matrix  elements  governing  ionisation  are  P  +A 

en  en ' 

between  bound  and  continuum  Born-Oppenheimer  eigen¬ 
functions  . 

In  chapter  I,  I  described  the  physically  unrealistic 
behaviour  of  the  uncorrected  pss  matrix  elements.  I  also 
discussed  how  Thorson  and  coworkers found,  for  H  +,  that 
it  is  possible  to  choose  a  switching  function  for  each 
bound  state,  such  that  the  size,  number,  and  range  of  the 


■ 


. 


couplings  from  that  state  to  the  continuum  are  all  greatly 
reduced.  In  this  thesis,  I  confirm  and  considerably  ex¬ 
tend  these  results,  and  also  make  them  much  more  precise. 

If  the  f's  so  obtained  can  be  called  "best"  choices,  in 
some  sense,  then  they  could  be  said  to  provide  an  "optimum" 

way  of  isolating  the  (so  far  undefined)  "deformation"  part 
of  P(R)  . 

Our  method  for  choosing  switching  functions  will  be 

described  at  length  in  chapter  V.  Both  the  studies 

reported  in  this  thesis  and  the  earlier  work  of  Thorson  et 
,  lb  .  , . 

al.  indicate  clearly  that  different  f's  are  needed  for 
different  bound  states;  in  chapter  V,  I  shall  also  consider 
how  to  modify  the  theory  to  take  account  of  this. 

F.  Rotating  Frame  Co-ordinates 

Up  to  this  point,  I  have  assumed  that  the  electronic 

states,  0n(r;R) ,  and  Hamiltonian,  h^ ,  are  expressed  in 

terms  of  electron  co-ordinates,  r,  in  a  reference  frame 

whose  axes  are  fixed  in  space?  in  particular,  the  gradient, 

(Vr)£,  is  evaluated  with  r  fixed  in  this  frame.  However, 

as  is  well-known,  the  molecular  electronic  states  and 

Born— Oppenhe imer  Hamiltonian  are  normally  described  in  a 

rotating  reference  frame,  whose  polar  axis  coincides  with 

R;  let  us  denote  the  electron  co-ordinates  in  this  frame 
,  ■+• . 

by  r  .  In  the  rotating  frame,  the  transformed  Hamiltonian, 


h0 '  ,  and  wavefunctions ,  <j>  1  ,  depend  upon  r '  and  R  =  |r|. 


' 


* 


* 


but  do  not  depend  on  the  orientation  of  the  vector,  R,  in 
space. 

The  transformation  linking  the  components,  (x',y', 
z'),  of  r',  to  the  components,  (x,y,z) ,  of  r,  is7a 


X  1 

cos 0  cos  $  cos 0  sin$  -sinO' 

r 

X 

y’ 

= 

-sin$  cos  $  o 

y 

z  ' 

^sinO  cos  $  sin0  sin$  cosG 

✓ 

iz. 

where  R  has  polar  co-ordinates,  (R,0,$).  Note  that  the 
^  _ax^s  coincides  with  the  line  of  nodes  and  is  therefore 
perpendicular  to  the  plane  containing  the  z-  and  z'-axes. 

Using  equation  (11-58) ,  we  can  express  the  com¬ 
ponents  of  the  gradient. 


($R)j  <Pn(r;R) 


in  terms  of  operators  in  the  rotating  frame.32  34 
that 


Noting 


V?;R)  =  V  (?’  ;R)  '  (11-59) 

a  straightforward  application  of  vector  calculus  yields 
the  result  that 


i  / 


-ifi  [3<f>n(r;R)/3R]j  =  -ifi  [  3$n  ’  ( r  ' ;  R)  /  3R]  j 


(II-60a) 


. 


. 

.«  4  •» 


(  "  ,  ,  .  •  'J  '  : 


-ifi  [3V?;8)/30]r  =  -  £  *n'(?*;R), 


(II-60b) 


ifl  [30n(r;R)/34>]j  =  {Lx,sin0  -  Lz,cos0}  4>n'(r';R); 

(II-60c) 

ys  <*.  /\ 

where  Lx, ,  Lyl,  Lz ,  are  the  components  of  the  electronic 
orbital  angular  momentum  in  the  rotating  frame.  The 
components  of  the  vector,  P(R) ,  are  then  given  by7a'32 


Pjn(R)  =  pin  z'  +  P<L  x'  +  p$  y'  / 
J11  jn  jn  jn  1  ' 


(11-61) 


A  A  A 


where  x',  y',  z'  are  the  unit  axis  vectors  in  the  rotating 
frame,  and 


,R 


P  jn  =  ~ifl  <<p  j  I  0/3R)j,  Un>, 


(II-62a) 


P0  = 

jn 


-  R 


-1 


<0  .  |  L  ,  I  (p  >  , 
3  1  y ' 1 rn  ' 


(II-62b) 


P jn  "  +  R_1  {<<f>j|Lx'IV  "  cot0  <^|£z.l*n>}. 

(II-62c) 

For  convenience,  I  have  dropped  the  primes  on  <J>  '  (r;R). 

The  components  of  A (R)  are  simply  defined  by  the  components 
of  the  vector,  s,  in  the  rotating  frame. 

In  the  classical  trajectory  description,  the  motion 
of  the  nuclei  is  confined  to  a  plane,  and,  without  loss  of 
generality,  we  can  take  this  as  the  (x',z')  plane.  More¬ 
over,  to  compute  the  matrix  elements,  P0  and  we 

Dn  jn' 


' 


. 


require  in  both  cases,  only  the  matrix  elements  of  the 

ladder  operators,  L± • .  Therefore,  the  behaviour  of  the 

angular  couplings  will  be  fully  covered  by  a  study  of  a 

single  angular  component,  which  we  take  to  be  (P0  +AX ' ) 

.  ...  jn  jn  * 

his  coupling  is  just  the  quantum  mechanical  Coriolis 
coupling,  seen  by  an  electron  in  the  rotating  frame. 

Asymptotically,  the  components  of  P(R)+A(R)  are  just 
[cf.  equation  (11-55) ] 

Pjn  +  Ajn  =  "ifl  <(J)j  |  4>n  >,  (II-63a) 

J  J  J 

and 


+  A' 


jn 


(II-63b) 
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CHAPTER  III 


BORN-OPPENHEIMER  WAVEFUNCTIONS 
A.  Introduction 

In  this  chapter,  I  shall  describe  the  computation  of 
the  eigenfunctions  for  two-centre,  fixed-nuclei,  one- 
electron  systems — the  Born-Oppenheimer  wave functions . 

These  states  are  solutions  of  the  adiabatic  Schrtfdinger 
equation , 

2  2  2 

{  (ti  /2mQ)  Vr  +  e  (zA/rA+zB/rB)  +  e]  T  (r ;  R)  =  0,  (III-l) 

where  e  =  e'(R)  -  Z^Z^e  /R  is  the  electronic  energy.  e 
may  be  either  positive  (the  electronic  continuum)  or 
negative  (discrete  bound  states) . 

It  is  well  known  that  the  Schrfldinger  equation  (1) 
is  separable  in  prolate  spheroidal  co-ordinates,  (£,n,<j>), 
where  (see  figure  3) 

^  —  ^A"*"rB^///^/  ^  ~~  $  =  azimuthal  angle 

( III-2 ) 

1  f  C  <  00  -1  <  n  <  +1  0  <  <p  <  2  7T ; 
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Figure  3.  Spheroidal  Co-ordinates 

S  =  (rA+rB)/R/  H  =  (rA-rB)/R,  tp 


and  and  rB  are  the  distances  of  the  electron  from  the 
nuclei  A  and  B,  respectively.  The  elect 


ron  has  mass,  m^ , 


and  we  use  the  Bohr  radius,  aQ,  as  the  unit  of  length;  the 

energy  is  in  Rydbergs  (1  Ry.  =  e2/2aQ  =  0.5  au) .  Y  can  be 
written 


y(r;R)  =  x(0  S(n)  <M<j>)  , 


(III-3) 


and  the  separated  equations  for  the  fact 


ors  are 


2  2  2 
d  $/d<f>  +  m  $  =  0 , 


(III-4 ) 


d/dri  {  ( 1— n  )dS/dn)  +  {pn±c2r)2+A-m2/ (l-r|2)  }S  =  0,  (III-5) 


d/dC  ((£  -l)dX/d£}  +  (q£-c2£2-A-m2/ ( £2-l)  }x  =  0;  (III-6) 


where  the  parameters,  p,  q,  and  c2,  are  defined 


p  =  R(Zb-Za) ,  q  =  R(Zn+Z J , 


B  A' 


(III-7a) 


2  .  2 
c  =  -  h  e  R  . 


(III-7b) 


A  and  m  are  the  separation  constants;  m  is  the  modulus  of 
the  azimuthal  quantum  number.  In  equations  (III-5)  to 
(HI  — 7)  ,  where  a  "±"  appears,  the  upper  sign  gives  the 


bound  state  equations,  while  the  lower  gives  those 


•  r ' 


* 


* 
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corresponding  to  the  continuum. 

In  the  united  atom  limit,  equation  (III-5)  tends  to 
the  equation  for  the  associated  Legendre  function,  P  m(ri)  , 
with  A  =  L(L+1);  and  equation  (III-6)  becomes  the  radial 
equation  for  either  the  spherical  Coulomb  function,  or  the 
associated  Laguerre  function.36  The  wavefunction  for  the 
azimuthal  co-ordinate  is,  of  course,  just 

$(<}>)  =  ( 2 7T )  exp(±inuf)),  with  m  =  0,  1...  (III-8) 

In  sections  B  and  C,  I  will  discuss  the  solutions  of 
equations  (HI-5)  and  (III-6)  ,  for  bound  and  continuum 
electronic  states,  respectively. 

B.  Bound  States 

1.  Background 

The  solution  of  equation  (III-l)  has  received  con¬ 
siderable  attention  for  negative  energies.  Subject  to  the 
usual  conditions  that  the  wavefunction  be  square— integrable , 
continuous,  and  have  continuous  first  derivatives,  the 
coupled  equations,  (III-5)  and  (III-6) ,  have  simultaneous 
solutions  only  for  discrete  pairs  of  values  of  the  par¬ 
ameters,  A  and  c.  In  the  limits,  R-*0  and  R->°°,  simple 
closed  formulae  exist  for  the  energies  and  wavefunctions , 
but,  in  general,  we  must  solve  the  equations  numerically. 

Various  authors  have  complied  extensive  tables  of 
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the  eigenparameters,  both  for  h/,37'40  and  for  unequal 
nuclear  charges.41,42  The  formal  properties  of  the  sol¬ 
utions  have  been  examined  in  great  detail  by  Power.43  To 

compute  solutions,  I  used  a  method  based  on  that  given  by 
Bates  and  Carson41  for  HeH++. 

For  a  given  bound  state,  the  calculation  of  the 
eigenvalues  of  A  and  c  proceeds  as  follows.  At  each 
consecutive  R-value, 

i.  Guess  approximate  values  for  the  eigenvalue,  c  ,  and 

g 

separation  constant.  A 

g 

ii *  Guess  a  spread,  6,  such  that  the  true  eigenvalue  is 
expected  to  lie  between  c  —  6  and  c  +5 

g  g 

iii.  For  a  series  of  values  of  c,  between  c  -6  and  c  +  6 

g  g  ' 

compute  the  value  of  A  which  gives  the  regular  sol¬ 
ution  to  equation  (III-5) ,  the  angular  equation 

iv.  Use  this  set  of  pairs  of  values  of  A  and  c  to 

determine  A  as  a  suitable  polynomial  function  of  c 

v.  Using  this  A(c) ,  compute  the  value  of  c  which  gives 
the  regular  solution  to  equation  (III-6) ,  the  radial 
equation.  This  is  the  desired  eigenvalue,  c 

vi.  Now  go  back  to  the  angular  equation,  and  compute  the 
separation  constant  corresponding  to  this  value  of  c 
Once  A  and  c  have  been  found,  the  corresponding  eigen¬ 
functions  can  be  computed  fairly  easily. 

In  the  next  two  sections,  I  shall  describe  the 


computation  of  solutions  to  the  angular  and  radial 


.  .  '  ■  F 


-  , 


equations  in  more  detail. 
2.  Angular  Equation 


d/dn  {(i-n2)ds/dn}  +  {Pn+c2n2+A-ra2/(i-n2) }  s  =  0.(111-9) 

We  require  solutions  to  (III-9)  which  are  regular  in 

the  interval,  -1  <  n  <  +1.  Following  Bates  and  Carson,41 
we  write 

s(n)  =  exp(-cn)  2  dL  PLm(n)  ,  (m-io) 

L=m 

where  (n)  is  the  L '  th  associated  Legendre  function. 
Substituting  (III-10)  into  (III-9) ,  and  using  the  recursion 
relations  satisfied  by  the  PLm,s,  we  obtain,  after  some 
manipulation , 

Z  dL  {PL-1  (P“2cL)  (L+m)/(2L+l)  +  PLm  [ c2+A-L  (L+l)  ] 

+  Pl+i  [p+2c  (L+l)  ]  (L-m+1) /(2L+1)  }  =  0.  (III-ll) 

Since  the  PLm's  form  a  complete,  orthogonal  set,  equation 

(III-ll)  can  only  be  satisfied  if  the  coefficient  of  P  m 

L 

vanishes  identically,  for  each  L.  This  leads  to  a  three- 

term  recursion  formula  linking  the  d  ' s: 

Li 


. 


; 


jt 
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dL-l  (P+2cL)  (L-m)  /  (2L-1)  +  dL  [c2+A-L  (L+l)  ] 

+  dL+l  fP“2c  (L+l)  ]  (L+m+1)  /  (2L+3)  =  0.  (III-12) 

For  very  large  L,  the  coefficients  satisfy  either 

dL+l  ~  "  (L/c)  dL  (HI-13a) 

or 


dL  -  +  (c/L)  dL-1 


(III -1 3b ) 


Clearly,  we  require  solutions  which  behave  like  (III-13b) ; 
for  a  given  c,  these  only  exist  for  particular  values  of  A 
-  the  eigenvalues.  To  determine  these,  I  used  Ince's  well 
known  method  for  finding  the  eigenvalues  of  a  tridiagonal 
matrix,44  as  described  by  Stratton  et  al.45 

For  computing  the  function,  S(n),  it  is  not  con¬ 
venient  to  use  the  expansion  given  in  (III-10) ;  convergence 
of  the  series  is  often  poor,  even  at  moderate  internuclear 
separations,  and  the  evaluation  of  integrals  over  n  is 

awkward.  I  therefore  chose  to  use  the  approach  introduced 
46 

by  Helfrich  and  expressed  the  wavefunction  as 


£ 

L=m 


m 


(n) 


s(n) 


(III-14 ) 


■ 
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. 
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ft 


This  gives  a  five-term 


i 


recursion  formula  for  the  d  ' s : 

L 


[L_2  c  (L-m)  (L-m-l)/[  (2L-3)  (2L-1)  ]  +  d^  p  (L-m)  /  (2L-1) 


+  dL  {c  [2L(L+l)-2m2-l]/[  (2L+3)  (2L-1)  ]  -  L(L+1)  +  A} 


+  dL+l  P(L+rn+l) /(2L+3)  +  dL+2  c2  (L+m+2)  (L+m+1)  /  [  (2L+5) 


( 2L+3)  3  =  0.  (III-15) 

Ihe  technique  for  finding  the  coefficients,  d  ,  given  the 
eigenvalue.  A,  is  described  by  Wilkinson^  and  the  series 
is  truncated  when  the  d^'s  become  sufficiently 
small.  Of  course,  the  wavefunction  cannot  be  computed 
until  the  eigenvalue,  c,  has  been  found. 

Let  us  now  turn,  then,  to  the  solution  cf  the  radial 
equation.  We  will  assume  that,  in  the  vicinity  of  the 
desired  eigenvalue,  A  is  known  as  a  function  of  c;  in  my 
programme,  I  chose  to  fit  the  computed  separation  constants 
to  a  cubic  polynomial  in  c,  by  the  least  squares  method. 

3.  Radial  Equation 

d/d£  { ( £2-l) dX/d£ }  +  {q£-c2£2-A-m2/ ( £2-l) }  x  =  0.(111-16) 


At  large  £,  the  two  linearly-independent  solutions 
to  equation  (III-16)  behave  like  exp(±c£).  Clearly,  we 
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are  interested  in  the  exponentially  damped  solution,  and, 

once  again  following  Bates  and  Carson,41  we  use  the  trans¬ 
formation  due  to  Jaffe: 

X(5)  =  U+D^c-l  { U-l) / U+l)  }m/2  exp ( — c 5 ) 

E  a  {  (£-l)/(£+l)  }n(III-17) 
n=0  n 

If  we  substitute  (HI-17)  into  (III-16),  collect  terms  in 

powers  of  (£-l)/(£+l),  we  obtain,  after  a  lot  of  algebra, 

a  three  term  recursion  formula  linkinq  the  a  '  s 

n  ' 

an-l  (n+m-q/2c) (n-q/2c)  +  aR  { ( 2n+m+l) (q/2c-2c-l) -2n (n+m) 

+q-A-c2}  +  an+1  (n+1)  (n+m+1)  =  0,  (III-18) 

with  a_x  =  0,  and  aQ  =  1. 

Given  A(c),  Ince's  method44  yields  the  desired 
eigenvalue,  c,  and  the  coefficients,  a^,  in  equation  (HI- 
17)  .  The  eigenenergy  is  then  given  simply  by 

e  =  -  (2c/R) 2 .  (III-19) 

4.  Normalisation 

The  wavefunction  is  normalised  in  the  usual  way. 


l  *  e 
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<V| y>  =  1. 


(III-20) 


The  volume  element  in  spheroidal  co-ordinates  is 

“  (R/2)  ( £  -r)  )  d£dr)d$,  (III-21) 

and  hence,  we  have  to  evaluate 

N2  =  (R/2)  3  {/^S2(n)dn  f°^X2(OZ2 d£  - 

f-is  (n)n2dn  /1x2(^)d^}.  (m-22) 

I  chose  to  normalise  the  coefficients,  d  ,  such  that 

Zj 

/_]_s  (n)dn  =  1.  (in-23) 

The  n-integrals  can  be  evaluated  analytically,  using  the 
orthogonality  relation  satisfied  by  the  Legendre  functions 

p-jm(n)P  m(n)dn  =  2/(2j  +  l)  (j+m)  l/(j-m)  i  6.  ;  (III-24) 

A  J  11  jn 

and  we  find  that 

y2  =  /^S2(n)n2dn  =  £  2  (L+m)  !/(L-m)  i  {d  2  [2L-1+2 (L2-m2)  ] 

L=m 

(III-25) 

+  2  dL+1dL_1  [L(L+1) -m(m+l)  ]  }/[  (2L-1)  (2L+1)  (2L+3)  ]  . 


' 


■ 


Thus,  equation  (III-22)  reduces  to 

N  =  ( R/2 )  /xx  (£)  (^2-y2)  d£.  (in-26) 

This  integral  has  to  be  done  numerically,  and  I  chose  to 
evaluate  it  by  Gaussian  quadrature.  Since  integrals  of 
this  general  form  appear  several  times,  I  will  describe 
the  method  used  for  evaluating  them  in  some  detail.48 
They  can  all  be  written 

1(a)  =  (R/2)  3  / 1  exp(-aC)  g(£)  d£.  (III-27) 

For  example,  in  equation  (III-26) ,  a  =  2c.  If  we  change 
the  variable  to  x  =  a(£-l),  we  find  that 

1(a)  =  (R/2)  e  a/a  /Q  e  x  g(l+x/a)  dx;  (III-28a) 

this  can  be  evaluated  by  Gauss-Laguerre  quadrature,  giving 

^  N 

1(a)  -  (R/2)  e  a/a  Z  w.  g ( £ . ) ,  (III-28b) 

i=l  1  x 


where 


€i  =  1  +  xi/a. 


( III-2  8c) 


The  integration  points , 


xi' 


and  weights ,  w. 


are  given  in 
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Abramowrtz  and  Stegun49  for  various  values  of  N,  and  I 

used  a  28-point  truncation  of  the  64-point  formula, 

although  this  was  probably  over-cautious.  Exploratory 

calculations  showed  that  this  formula  will  produce  about 

12  (!)  figure  accuracy  for  a  wide  variety  of  functions, 

g(x)  ,  whereas  the  error  in  the  wavefunctions  is  about  one 
part  in  106. 

The  programme  computes  the  value  of  N,  and  stores  it 

for  later  use.  Let  us  now  turn  to  the  problem  of  computing 
the  continuum  wavefunctions. 

C.  Continuum  States50 

1.  Background 

The  electronic  continuum  states  for  the  two-centre, 

one-electron  system  have  received  very  little  attention, 

compared  with  that  given  to  the  bound-state  functions. 

Bates,  Opik  and  Poots  ^  describe  the  solution  for  H  +  ,  in 

connection  with  photo-ionisation  cross  sections  for  the 

molecule-ion.  Recently,  Ponomarev  and  Somov52  have 

described  the  solution  for  HeH++,  and  an  incorrect  solution 

corrected)  for  arbitrary  charges  has  been  published 

by  Greenland.  I  have  used  a  method  similar  to  that 

52 

of  Ponomarev  and  Somov.  Regular  solutions  exist  for  all 
positive  energies,  and  for  discrete  values  of  the  separ¬ 
ation  constant,  A. 


. 
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2.  Angular  Equation 


d/dn  ( (i-n2)ds/dn>  +  {pn-c2n2+A-m2/(i-n2) }  s  =  0.(111-29) 

We  require  solutions  to  equation  (III-29)  which  are 
regular  in  the  interval,  -1  <  n  <  +1.  Now,  (III-29)  is 
obtained  from  the  corresponding  bound  state  equation  if  we 
replace  c  by  ic.  This  suggests  that  we  can  solve  (III-29) 
in  a  manner  analogous  to  that  used  for  (III-9) ;  we  write 

S(n)  =  exp(-icn)  £  dL  PLm(n)  ,  (III-30) 

L=m 

and  we  obtain  a  three-term  recursion  formula  for  the  d  's: 

L 

dL-l  (p+2icL)  (L-m)/(2L-l)  +  (A-c2-L(L+1)  }  + 

(b+1)  }  ( L+m+ 1)  /  ( 2L+ 3)  =0,  (III-31) 

with  dm_!  =  °-  The  eigenvalues  of  A  are  real  (as  Ponomarev 
52 

and  Somov  point  out,  they  depend  only  on  the  absolute 
values  of  the  complex  coefficients  in  (III-31) ) ,  and  they 
can  be  determined  by  Ince's  method44  for  a  given  continuum 
energy  and  R- value. 

As  with  the  bound  state  equation,  the  expansion, 
(II1~30) ,  does  not  provide  a  convenient  representation  of 
the  wave  function ,  and  we  express  the  solution  as 


1 


: 


■ 


♦ 


7 


^  "  Jm^  h>»'  (HI-32) 

which  yields  a  five-term  recursion  formula  for  the  d  's: 

Xi 

dL_2  c2  (L-m)  (L-m-1)  /  [  (2L-3)  (2L-1)  ]  -  d^  p  (L-m)  /  (2L-1)  + 
dL  (L(L+1)  -A+c2[2L(L+1)  -2m2-l]/[  (2L-1)  (2L+3)  ]  }  - 

^0 

dL+l  P  (L+m+1)  /  ( 2L+3)  +  dL+2  c2  (L+m+2)  (L+m+1)  /  [  (2L+3)  (2L+5)  ] 

=  °*  (HI-33) 

Given  an  eigenvalue.  A,  the  dL's  can  again  be  computed  by 
the  method  described  in  Wilkinson.47  Notice  that  equations 
(III-29)  to  (III-33)  can  all  be  obtained  from  the  corres¬ 
ponding  bound  state  equations,  (111-9,10,12,14,15),  by 
replacing  c  with  ic. 

3.  Radial  Equation 

d/d£  {(£  -1) dX/d£}  +  (q£+c2£2-A-m2/ ( £2-l) }  x  =  0.(111-34) 

We  require  the  solution  of  (III-34)  which  is  regular 
at  ^  =  !•  At  large  £ /  this  solution  has  the  form 

req 

x  y(£)  ^  (B/r)  sin  {kr+yln  (kr)  +  6  }  , 


(III-35a) 


'  '  •  J  •  1 


'  ■  ■  ■' 1  ; :  >  J : ! 


where 


Y  =  q/  (2c)  =  (zA+zB)/k, 

cC  ^  kr,  k2  =  £. 

B  is  the  normalisation  constant,  and  5  is  the  phase  shift. 

Equation  (HI-34)  is  solved  by  starting  the  solution 
near  the  origin  with  a  series  expansion,  and  continuing 
this  solution  outward  by  numerical  integration,  until  it 
can  be  matched  with  asymptotic  solutions  to  determine  the 
phase  shift.  The  solution  regular  at  5  =  1  has  the  form 

X(5)  =  {(C-D/U+l)  }m/2  F(C),  (III-36) 

where  F(4)  satisfies  the  equation 

(4'  -1)F  '(4)  +  2  (m+4)  f  '  (4)  +  (c242+q4-A)  f  (4)  =  0.(111-37) 

Near  the  origin,  this  can  be  written  as  a  power  series  in 
(4-D  , 


(III-35b) 

(III-35c) 


Freg(4)  =  Z  a  ( 4-1)  n , 
n=0  n 


(III-38) 


and  this  yields  a  four-te 

coefficients,  a  : 

n 


rm  recursion  relation  for  the 


"I 


2  2 

an- 2  c  +  an-l  (2c  +c2)  +  an  [n  (n+1)  +c+q-A]  + 

2an+l  ^n+1)  (n+m+1)  =  0,  (III-39) 

where  a_2  =  a_ ^  =  0;  and  aQ  =  1. 

The  solution,  F  g(£),  is  continued  numerically  by 
integration  of  (TH-37)  using  a  variable-step,  Adams - 
Moulton-Bash forth,  predictor-corrector  method,54  until  a 
value  of  £  suitable  for  matching  Freg(£)  with  asymptotic 
solutions  is  reached.  However,  for  evaluating  integrals 
numerically,  we  have  seen  that  the  value  of  the  integrand 
is  required  at  specific  values,  ,  whereas  the  A-M-B 
algorithm  gives  the  function  at  successive  points,  separ¬ 
ated  by  a  predetermined  step-size.  The  programme  must 
thus  include  a  modification  to  evaluate  Freg(£.):  whenever 
the  A-M-B  integration  passes  one  of  the  ,  the  programme 
pauses  to  use  the  fourth-order  Runge-Kutta  method54  to 
calculate  and  store  the  Freg(£.). 

At  large  £,  asymptotic  solutions  to  (III-37)  are 

=  (£+1)  exp±i  {c£+yln  (£+1)  }  u“(£),  (III-40) 

and  if  we  write 

u±U)  =  1  b  1  {2/  [c  ( £+1)  ]}n, 

n=0 


(III-41) 


. 


■ 

- 

,  ‘a  ■  -  !  o ■  u  ' 


the  coefficients 


satisfy  a  three— term 


recursion  formula: 


bn-l  c{n(n+m)-y2  +  iy  (2n+m)  }  -  {n  (n+1) -Y2+q+c2-A 

+  iYt(2n+l)-2c(2n+m+l)]}  ±  4i  b^  (n+1)  ,  0,  (III.42) 

with  b _l  =  0,  =  1.  Evidently,  u'f?)  =  {u+(C)}*.  The 

expression,  (III-41),  is  an  asymptotic  expansion;  if  g  is 
large  enough,  it  yields  sufficiently  accurate  solutions, 

F  (5).  The  phase  shift  is  determined  by  matching  FregU) 
and  dFreg(U/dS  to  some  suitable  linear  combination  of  the 
asymptotic  solutions;  that  is,  we  write 

FrSg(C)  =  D  {eia  F+(5)  +  e-ia  F-(C)} 

(III-4  3a) 

CFreg(5)]'  =  D  {eia  [F+(C)]-  +  e-ia  rF~  (C)  )  •  > 
and  determine  a.  We  find  that 


tana  =  Re (h*) /Im (h1) 


(I II -4  3b) 


h+  =  (u±(C)  ' 


u  tw(C)  +  ig  (0  '  ]  >  exp±ig(£),  (m-43c) 


where 


w(C)  =  [Freg(C)  ] '/Freg(C)  +  (C+l)’1 


(III-44a) 


' 


r*  ■ 


■ 
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and 


g(S)  =  c£  +  yin  ( £+1)  . 


(III-44b) 


From  equations  (IIX-40)  and  (III-43),  it  follows  that 


,reg 


x  (5)  ^  fc/fs+l)}  sin {c£+yln  [ c  ( 5+1)  ]  +6  }  ,  (III-45a) 


6  =  a  -  yln(c)  +  n/2,  Y  =  q/(2c)  . 


(III-45b) 


In  the  united  atom  limit,  6  tends  to  the  appropriate 
Coulomb  scattering  phase, 


lim 

R+0 


6  =  aL  "  (Ltt/2 )  +  yin  (2)  , 


(III-4  6a) 


aL  =  Arg{r (L+l-iy) }. 


(III-4  6b) 


lb 

Levy  and  Thorson  have  reported  phase  shifts  for 
the  H2  problem,  but  these  contain  an  arithmetic  error. 
The  correct  phases  are  obtained  by  subtracting  c  =  J$kR 
from  their  values;  other  quantities,  such  as  the  matrix 
elements  they  reported,  are  unaffected.  This  error  has 
also  been  noted  by  Greenland  and  Greiner. 53b 

4.  Normalisation 

The  wavefunction  is  normalised  so  that  the  density 
of  continuum  levels  is  included. 


. 
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^e'L’m'  1 4'eLm>  “  6  U  ' _e)  6L.  L< « 

■Li  ij  m  m 


(III-47) 


If  the  coefficients  in  the  angular  eigenfunction 
normalised  such  that 


are 


s2 (n )  dn 


this  requires 


=  1, 


(III-48) 


C  2/(2ttcR)^  -  2/[R(irk)  . 


(III-49) 


If  one  wishes  to  express  the  energy  in  atomic  units  (2  Ry. 

=  1  au) ,  one  has  to  be  rather  careful  at  this  point.  The 
density  of  states,  dn/ds,  in  atomic  units  is  twice  that  in 

y  ergs,  and  hence,  Cau  =  2 **  .  The  programme  computes 

c  (m  Rydbergs)  and  stores  it  for  later  use. 

5.  Remarks 

The  representation  of  the  solutions  differs  in  some 
ways  from  that  given  by  Ponomarev  and  Somov. 

a.  I  use  an  expansion  in  Legendre  polynomials  to  compute 
the  angular  solution,  whereas  they  employ  a  power 
series  near  £  =  -1,  followed  by  numerical  integration. 

b.  Using  the  extended  asymptotic  solutions  for  the  radial 
wavefunction  (equations  (III-40)  to  (III-42))  allows 
more  rapid  calculation  of  both  phase  shifts  and  wave- 
functions.  The  point,  £*,  at  which  the  numerical 


«  • 


■ 


solution  is  matched  to  the  asymptotic  solution,  is 
given  by 


★ 

~  3Y;  (III-50a) 

Ponomarev  and  Somov52  state  that  they  found 

c£  -  50  [a-c2+y2|.  (III-50b) 


c.  The  representation  of  XregU)  near  C  =  1  is  also 

slightly  different,  but  this  is  not  significant. 

53a 

P.  T.  Greenland  proposes  a  solution  for  the 
continuum  states  which  is  not  correct.  In  particular: 

a.  He  imposes  constraints  which  imply  that  the  states, 
S(n),  are  parity  eigenstates,  regardless  of  the  charge 
asymmetry;  this  is  clearly  false 

b.  He  attempts  to  compute  Xreg(£)  by  fixing  the  asymptotic 

behaviour  of  the  solution,  rather  than  the  behaviour 
at  £  =  1 


Let  us  now  look  at  some  representative  results. 

6.  Results 

The  notation  we  use  labels  the  continuum  states  with 
the  energy,  c,  the  united-atom-limit,  angular  momentum 
quantum  number,  L,  and  m,  the  modulus  of  the  azimuthal 
quantum  number. 

Table  1  gives  some  representative  values  for  the 
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separation  constants  and  phase  shifts  in  HeH++;  data  shown 
are  for  e  -  1.0  Ry. ,  and  the  states,  so  (L=m=0)  ,  pa  (L=l, 
m=0),  and  Ptt  (L=m=l)  ,  for  internuclear  separations,  0<R<12 
au.  In  the  united  atom  limit,  the  phase  shifts  all  agree 
with  the  appropriate  Coulonfc  phases  ((III-45)  and  (HI-46)) 

modulo  2tt.  Tables  2  to  4  illustrate  the  behaviour  of  the 

expansion  coefficients  H  , 

'  dL'  t^e  an<?ular  eigenfunctions 

for  the  same  states,  for  0<r<6  au. 

Figure  4  shows  the  separation  constants  for  the  model 
system  with  (Zft/ZB)  =  1.2,  at  energy  e  =  0.5  Ry.,  for  the 

states,  so,  pa,  pr,  da,  dr,  d<5 .  The  results  may  be  applied 
directly  to  a  real  bare-nucleus  system  with  z  ' ,  z  ■ 

integers  (e.g.,  12,  10),  by  scaling  all  distances  as 

(ZB')  ,  all  energies  as  (Z  ' )  ^ . 

£ 


Figures  5  and  6  depict  phase  shifts  for  the  model 
system  at  e  =  0.5  Ry. ,  and  for  HeH++  at  e  =  1.0  Ry. ,  for 
the  states,  pir,  da,  dr,  fa,  d6 ,  fr ,  ga ,  fS,  gir,  ha.  It  is 
interesting  to  note  the  systematic  interleaving  and 
ordering  of  the  phases  in  the  large  R  region.  They  appear 
to  be  grouped  according  to  the  value  of  (L+ra) ,  with  phase 
shifts  decreasing  as  that  index  increases;  within  any  group 
of  given  (L+m) ,  they  are  ordered  by  the  value  of  L,  with 
phase  shifts  decreasing  as  L  increases.  Members  of  a  group 
never  cross  each  other.  This  behaviour  seems  to  be  quite 
general  for  all  charge  ratios  and  energies. 
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TABLE  1 

SEPARATION  CONSTANTS  AND  PHASE  SHIFTS  FOR  HeH++; 


e  =  1.0 

Ry.  and  6  IS 

IN  UNITS  OF 

7 r 

sa 

WAVE 

pa 

WAVE 

P7T 

WAVE 

R 

A 

6 

A 

6 

A 

6 

0.0 

0.000000 

0.32661 

2.000000 

1.42903 

2.000000 

1.42903 

0.5 

-0.019726 

0.21542 

2.061212 

1.48671 

2.000083 

1. 39992 

1.0 

-0.067583 

0.08905 

2.231426 

1.54631 

2.001287 

1. 33885 

1.5 

-0. 119657 

-0.02076 

2.480920 

1.52845 

2.006210 

1.26945 

2.0 

-0. 154305 

-0.11892 

2.780038 

1. 47423 

2.018430 

1.19966 

2.5 

-0.157431 

-0.20910 

3. 106771 

1.40952 

2.041735 

1.13187 

3.0 

-0.121627 

-0.29334 

3. 448244 

1. 34357 

2.079508 

1.06667 

3.5 

-0.044246 

-0.37280 

3.800033 

1.27934 

2.134352 

1.00411 

4.0 

0.074055 

-0.44815 

4. 164259 

1.21740 

2.207954 

0.94408 

4.5 

0.230295 

-0.51974 

4.547059 

1.15752 

2 .301088 

0.88642 

5.0 

0.419874 

-0.58781 

4.956207 

1.09926 

2.413727 

0. 83100 

5.5 

0.637180 

-0.65248 

5.399372 

1.04218 

2.545180 

0 . 77770 

6.0 

0.876168 

-0 . 71384 

5.882991 

0 .98597 

2.694259 

0.72643 

6.5 

1. 130917 

-0.77200 

6.411565 

0.93045 

2.859430 

0.67709 

7.0 

1. 396101 

-0 . 82707 

6.987263 

0.87558 

3.038956 

0.62962 

7.5 

1.667318 

-0.87921 

7.609767 

0. 82142 

3.231025 

0.58394 

8.0 

1.941226 

-0 .92858 

8.276389 

0.76813 

3. 433846 

0.53999 

8.  5 

2.215513 

-0.97539 

8.982417 

0.71591 

3. 64  5729 

0.49768 

9.0 

2.488734 

-1.01984 

9. 721666 

0.66497 

3.865132 

0.45694 

9.5 

2.760112 

-1.06213 

10.487144 

0.61552 

4.090694 

0.41771 

10.0 

3.029326 

-1. 10245 

11.271748 

0.56773 

4.321242 

0. 37990 

10.5 

3.296346 

-1. 14098 

12.068871 

0.52174 

4 . 555787 

0.34345 

11.0 

3. 561303 

-1.17788 

12.872869 

0.47758 

4 . 793513 

0.30828 

11.5 

3.824407 

-1.21327 

13.679304 

0.43529 

5.033756 

0.274  32 

12.0 

4.085888 

-1.24730 

14.484988 

0. 39480 

5. 275984 

0.24149 

I 
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000817  .000627  -.000054  -.000015  .000002 


TABLE  3 

COEFFICIENTS  OF  THE  n-EXPANSION  (EQUATION  (111-32)  OF  TEXT) 
FOR  HeH++,  pa  WAVE,  AND  e  =  1.0  Rv 


82 


.. 


4 


i 

83 


CM 


Figure  4.  Continuum  separation  constants,  for  s,  p, 


and  d  waves 


R  (a.u.) 


3rr/2 
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Figure  5.  Phase  shifts  for  model 


system 


-  3tt/2 


3  tt/2 
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Figure  6.  Phase  shifts  for  HeH++ 


-3n/2 


CHAPTER  IV 


COMPUTATION  OF  MATRIX  ELEMENTS 

A.  Introduction 

In  chapter  II ,  we  saw  how  to  formulate  the  theory  of 
slow  atomic  collisions  in  a  way  which  takes  account  of 
electron  translation.  I  showed  that  the  nonadiabatic 
couplings  which  appear  in  pss  theory  must  be  modified  by 
important  corrections.  These  corrections  exactly  cancel 
the  spurious  asymptotic  pss  couplings,  and  also  produce 
substantial  reductions  in  the  effective  coupling  matrix 
elements  at  finite  internuclear  separations.  The  deriv¬ 
ation  is  approximate,  since  it  assumes  that  the  heavy 
particles  behave  classically,  but  the  main  conclusions  are 
the  same  as  those  in  the  rigorous,  quantum  mechanical 
derivation,  given  by  Thorson  and  Delos. 

The  formulation  employs  a  switching  function,  f(r;R), 
to  represent  the  translational  motion  of  the  electron  as  a 
local  function  of  its  position.  Apart  from  the  asymptotic 
constraints  (f  tends  to  +1  on  B,  and  -1  on  A) ,  f  is  un¬ 
defined,  and  the  derivation  provides  no  criteria  for  its 
specification.  The  problem  of  choosing  f  remains  unsolved. 
Many  of  the  obvious  defects  of  pss  theory  are  removed  by 
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any  reasonable  choice  of  f  which  meets  the  asymptotic  con¬ 
straints,  but  there  are  many  transitions  for  which  the 
corrected  matrix  elements  are  very  sensitive  to  that 
choice.  Studies  of  the  f-dependence  can  provide  useful 
insight,  particularly  if  the  behaviour  of  the  matrix 
elements  follows  some  obvious  pattern. 

The  corrected  couplings  for  transitions  from  low- 
lying  bound  states  to  the  electronic  continuum  are  known2 
to  be  highly  sensitive  to  variations  in  f.  In  this  chap¬ 
ter,  I  will  describe  the  computation  of  matrix  elements 
for  these  transitions.  The  motivation  for,  and  the  results 
of,  these  calculations  will  be  discussed  in  detail  in  the 
next  chapter. 


B.  Coupling  Matrix  Elements 

In  pss  theory,  the  matrix  elements  for  nonadiabatic 
coupling  are  defined 

Pjn(R)  =  <(^  j  I  (^R)  ?  I  c^n>  '  (IV-1) 

where  <£n(r;R)  are  the  Born-Oppenheimer ,  molecular  elec¬ 
tronic  states;  R  is  the  internuclear  co-ordinate,  and  r  is 
the  electron  co-ordinate,  measured  from  the  centre  of  mass 
of  the  nuclei  (CMN) .  P^n  has  "angular"  as  well  as  "radial" 
components,  since  the  molecular  states  are  defined  in  a 
rotating  frame,  in  which  Coriolis  forces  appear. 


. 
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In  chapter  II,  I  showed  that,  in  the  corrected 
theory,  P(R)  is  simply  replaced  by  Pg(R)  +  Ag(R),  where 


and 


(im/h)  «t'j|[hB0,sg]|<t,n>, 


(IV-3) 


with 


=  h  f(r;R)  r  ; 

9 


(IV-4 ) 


hB0  is  the  Born-Oppenheimer  electronic  Hamiltonian,  and  r 

g 

is  the  electron  co-ordinate,  measured  from  the  geometric 
centre ,  rather  than  the  CMN. 

is  sufficient  to  compute  the  radial  component  and 
one  angular  component  of  the  above  couplings.  The  explicit 
expressions  used  are  as  follows, 
i .  Radial  Couplings 


tE.9]jn  ”  <(t>j  I  (  9Ve/9R)  +  Un>/(en-ej),  (IV-5a) 

[^9]jn  =  (im/2ti)  (e_.-en)  ^  j  I  zgf  Un>  '  (IV-5b) 


where  en  are  the  Born-Oppenheimer  electronic  eigenvalues, 


. 
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Zg  is  the  component  of  r  on  the  molecular  axis,  and  V 

^  e 

the  electronic  potential  energy.  To  obtain  (IV-5a)  from 
(II-62a) ,  I  have  used  the  Hellmann-Feynman  relation. 

We  need  to  express  OVe/9R)j  in  terms  of  prolate 
spheroidal  co-ordinates  (see  figurl  3).  We  find  that 


(  3V  / 3R)  -*• 
e  r 


=  (e“/2) 


{ZAcosVrA' 


-  zBcos0B/rB  }.  (IV-6) 


But 


COS0J  -  (5n-fJ)/(5-fJn),  (iv-7) 

where  J  is  A  or  B,  and  fft  =  -1,  ffi  =  +1.  Thus,, we  find 
(R/2)3(52-n2)  ove/3R)j;  =  (Re2/ 4)  {z  [n  ( s2-l)  +5 ( 1-n2)  ] / 

g 

(C+n)  -  ZB t n  ( ?2-l) ( i— n 2 )  ]/( ?-n) 2 } .  (iv-8) 


ii c  Angular  Couplings 


tP9]jn  =  "  R_1  <*jlLy9l*n>'  (IV-9a) 

[Ag)®n  =  (im/2fi)  (£j-£n)  <*j|x  f|*n>,  (IV-9b) 


where 

ecular 


x  is 

g 

axis 


the  component  of  perpendicular  to  the  mol- 

A 

and  L  ^  is  the  y-component  of  the  electronic 


■ 


,  0%j  i:  '-\  •' 


« 


orbital  angular  momentum,  with  respect  to  the  geometric 
centre.  Now 


Lyg  ~  “  (i/2)  <L+9  "  L_9)  /  (IV- 10a) 

A 

where  L±g  are  the  ladder  operators  raising  and  lowering 
the  molecular-axis  component  of  orbital  angular  momentum; 


L±g  =  &  exp  (±i$)  {±8/80  +  i  Cot0  8/80}.  (IV-lOb) 


Using  the  relations, 


rgCos 0  -  (R/2 Hn;  rgsin0  =  (R/2)  (^2- D^d-n2)35, 
it  can  be  shown  that 

A 

L±9  =  fi  expf  +  ii}.)  {±[(5  -l)’5(l-n2)*5(r|3/35-53/3n)/(C2-n2)  ] 
+  i  5n/[  (52-l)^(l-n2)!s]  3/3<j) } .  (iv-ii) 


Computation  of  matrix  elements  reduces  to  separate 
quadratures  over  n,  the  "angle"  variable,  and  £,  the 
"radial"  variable.  The  n-integrals  were  performed  either 
analytically ,  or  by  32— point  Gauss— Legendre  quadrature, 
whichever  was  more  convenient;  and  the  ^-integrals  were 


done  by  a  28-point  truncation  of  the  64-point  Gauss- 


. 


■  : 


.  .  . 


. 


* 


Laguerre  quadrature  formula.  All  integrals  are  straight¬ 
forward,  except  for  those  of  (3V /3R)-  ,  where  a  logarith- 

G  IT 

g 

mic  singularity  must  be  extracted  analytically.  I  will 
describe  the  details  of  the  matrix  element  computations  in 
the  next  three  sections.  Unless  the  calculation  of  integ¬ 
rals  holds  a  particular  fascination  for  the  reader,  I 

suggest  that  he  pass  on  quickly  to  the  summary  given  in 
section  F. 

For  those  still  with  me,  I  will  first  establish  the 
notation.  To  avoid  confusing  initial  and  final  electronic 
states,  this  is  slightly  different  from  that  used  in 
chapter  III. 

i.  Bound  States 


IV  -  srm(rl)  W?)  '  (IV-12) 

where  n  and  r  are  the  united-atom-limit,  principal  and 
angular  momentum  quantum  numbers,  respectively,  and  m  is 
the  modulus  of  the  azimuthal  quantum  number. 


Srm(rt)  =  Z  dLrm  PLm(rl)' 

L=m 

Xnrm(5)  =  exP('cn?)  Fnrm(5)' 


c 

n 


(IV-13a) 


(IV-13b) 


h  e 


n 


where  c^  is  the  eigenparameter , 


2 


* 


■ 


'  '  H 


"  V  V. 
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ii.  Continuum  States 


<(^j  I  Xe  '  r  '  m '  ^  ^m'  *  (IV-14) 

T  and  m  have  the  same  meaning  as  the  corresponding 

unprimed  parameters,  and  s'  is  the  continuum  electronic 
energy . 


S 


r 


xre9 

eT'm' 


(O 


Z  dL 
L=m '  L 


r'm' 


PLm  <n>. 


=  t(5-l)/(S+l)]m/2  (5) 


(IV-15a) 

(IV-15b) 


I  will  assume  that  both  the  bound  and  continuum 
states  have  been  appropriately  normalised. 


C.  Computation  of  the  Matrix 

Elements  of  L 

— — - y 


1.  Introduction 

A 

The  operators,  L±g,  link  states  with  m'  =  m±l,  and 
the  integration  over  <f>  can  be  done  by  inspection.  Thus, 
in  effect,  the  operator  is 


L±g  =  fi  (±8/30  -  m  cot0).  (IV-16) 


All  the  integrals  over  n  can  be  performed 
with  the  help  of  the  recursion  relations 
PLm(n) 's,  and  integration  by  parts.  When 


analytically , 
satisfied  by  the 
<4>j  I  is  a 


i  ■ 


- 


..  -  ■  s\  ;  .it  “ 


- 


94 


continuum  state,  the  ^-integrals  are  all  of  the  form 

I(cn)  =  *  1  exP(“cnC)  g  ( 0  d£  (IV-17a) 

-  exp(-cn)/cn  Z.  w.  g(C.)  (IV-17b) 

h  =  1  +  W 

where  and  are  the  Gauss-Laguerre  quadrature  points 
and  weights. 

I  will  give  the  results  for  initial  states  with  m  = 
0  (a  states)  ,  and  m  =  1  (tt  states)  . 

2*  Ig_>Tr  =  <e  '  F  '  tt  |  L+g  |  nTa>  (m=0,  m'=l) 

From  equation  (IV-11) , 

■^a+Tr  =  ^  ^R//2^  f ]_  xe; Tp  .^(O  (a(a-*-Tr)  d/d£ 

-  B(a- tt)  C>  Xnra(0  d£,  (IV-18a) 


where 


a(a-*7r)  -  s r ' tt ( n )  d-n2)35  n  sra(n)  dn  (iv-i8b) 

=  (6/5)  d?T  71  d-.  +  Z  2(L+1)/[(2L+1)  (2L+3)  ] 

L=1 

^dLr  11  a  To  L/  ( 2 L—  1 )  +  (L+2)  /  ( 2L+5)  [  (L+3)  d-^'^d  Fa 

Lt  2.  L 
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-  L  d. 


tt,  TalX 
dL+2  ^  ^ ; 


and 


B(a^7T)  =  Srl7T(n)  (l-n2)^  dSrcr(n)/dn  dn  (IV-18c) 

=  Z  dTr,,rdr  ra  2L(L+1)/(2L+1)  . 

L=1  L  L 

Hence , 

(R/2)3  exp  (-cn5)  X^.^C)  (g2-!)1* 

dFnra(^)/d^  ”  (acn+6C)  F^U)}  dS-  (IV- 19) 

3.  l7T_).a  =  <e  '  T  '  a  |  L_^ |  nr7T>  (m=l ,  m'  =  0) 

If  we  notice  that 

Cn  ( S2-n 2 ) / [  (C2— D^d-n2)15]  =  S(S2-l)  Vd-n2)*5  + 

nd-n2)!s5/(52-l),s, 

then,  for  m  =  1,  after  integration  over  <p , 

U2- n2)  L±g  =  ti  {[±(C2-I)!s3/3C  -  ?/(?2-l)  **  ]  n  d-n2)  **  - 
(±d-ri2)!53/3n  +  n/(l-n2)!s]5(52-D!s}. 


(IV-20) 


■ 


*  ' 
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After  some  manipulation,  we  obtain 

Iir-*-a  =  "  R  (R/2’3  '1  X^^,a(C)  {a  ( TT^-a )  d/d$  (i^-l)*5  + 

B(ir-a)  C  (C2-l)  X^U)  d$,  (IV-21a) 

where  (cf.  equation  (IV-18b) ) 

a(w-a)  =  (6/5)  d0r’ad2r,r  +  Z  2  (L+l) /[  (2L+1)  (2L+3)  ] 

L=  1 

{dLr  l/(2L-1  )  +  (L+2)/(2L+5)  [  (L+3)  dTr'adTr!!  - 

L  L+  2 

L  dL+2adLrTT^;  (IV-2 lb) 

and 


B(TT^a)  =  Z  dLr'adLr7r  2L(L+1)/(2L+1)  .  (IV-21c) 

L~  1 


Hence , 


I 

TT~+0 


ft  (R/2)  3  /“  exp(-cn5)  X^f^,a(5)  {6C-acn+ad/d5} 


(^-1)^  FnF7r(5)  d£.  (IV-22) 

4*  =  <e  T  '  6  |  L+g|  nT7r>  (m=l,  m'=2) 

From  equation  (IV-20) ,  we  see  that 


' 


.  •  : 


1  1  ' 
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^  (R/2)  ^ a  (it-+6 )  [  ( £2-l)  **  d/d£  - 

e/tc2-!)*5]  -  6  (tt-i-6  )  £  (£2-l)  Xnr7T  ( c)  d£,  ( IV-2  3a) 

where 

a(TT-6)  =  /+]■  srlfi(n)  n  (l-n2)*5  sFir(n)  dn 

=  (16/7)  d,r  <SdiriT  +  £  2L(L+1)  (L+2)/[  (2L+1)  (2L+3)  1 

L=2 

{3dLF  (L-1)/(2L-1)  +  (L+3)/(2L+5)  [  (L+4)  d/'^d  r,T  - 

Li+  2  L 

(L-i)  dLr  5dL^]};  ( IV-2 3b) 

and 

6(tt^6)  =  f+i  sr,6(n)  (( l-n 2)i*  d/dn  +  n/d-n2)^}  srTT(n)  dn 

r  '  6  Ftt 

=  f  dL  dL  2(L+2)  (L+l)  L  (L-l)  /  (2L+1)  .  (IV-23c) 

L — 2 

There  fore  , 

Iwh-6  =  (R/2)3  exp(-cnC)  X^,6(5)  {(C2-l)!s[ad/dC-acn- 

60  -  aS/U2-!)15}  Fnr7[(5)  d£. 


(IV-24) 


■ 


* 


- 


•• 

•  i ' 
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This  completes  the  calculation  of  the  matrix  elements 

A 

of  L±  .  The  angular  component  of  P(R)  is  then  given  by 


g  1 0 


jn 


±  (i/2R) 


+ 

■^m+m '  * 


(IV-25) 


We  now  turn  to  the  matrix  elements  of  ( 9V  /3R)-> 

e  r  * 

g 

D.  Matrix  Elements  of  0/9R)  + 

— - — - - — i - r 

g 

1.  Introduction 
We  have  to  evaluate 


Jm-»-m '  =  <e ' r 'm' I (  3Ve/9R)  j  |nrm>, 

g 


(IV-26) 


where  (9Ve/9R)+  is  given  by  equation  (IV-8) .  The  integral 

t  g 

over  (p  is  trivial  and  yields  the  selection  rule  that  the 
matrix  elements  are  non-zero  if,  and  only  if,  m'  =  m. 
Therefore , 


!m  =  (Re  / 4 )  /“  Xnrm(C)  {ZA  I+(?)  -  Zfi  X_(?)  }  d5, 


(IV-2  7a) 


where 


i±(5)  =  ±  U(l-n2);n(C2-l)}/(C±n)2  sr,m(n)  srm(n)  dn. 

(IV-2  7b) 

If  we  factor  Sp,^(ri)  and  Sp^(r|)  into  even  and  odd  parts, 

s  (n)  =  se(n)  +  s°(n)  , 


. 


- 


■ 


.  rs 
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we  find  that 


V«  =  1  f5(l-n2)  -n  (C2-i) }/ ( C-n ) 2  {s* m  } 


{Srm(n)+Srm(n)  }  dn-  <IV_28) 


Hence , 

im'(5)  =  R  {zA  I+(C)  -  zB  i_ (5)  } 

plSf'nS?n,*Sr'niSfn,1  1  dn’  dV-29a) 

where 

q  =  R(ZB+ZA) ,  p  =  R(ZB-ZA).  ( IV-29b) 

Thus ,  we  have  to  compute 

Im  =  (e2/4)  fl  Xe^r'm(?)  Xnrm(C)  V  (5)  d?-  (IV-29c) 

The  integrals,  (IV-29),  are  not  trivial,  and  it  is 
instructive  to  examine  their  general  behaviour,  by  looking 
at  simpler  integrals  which  have  similar  properties. 


. 


f  « 

.. 


2.  Model  Problem 
Consider  the  integral 


100 


i(C)  =  {^(l-n  )  -n  (c2-i)  }/ ( C— n )  2  pL(n)  dn.  (iv-30) 

If  we  use  the  fact  that 

Qn(x)  =  Pn(y)/(x-y)  dy,  (IV-31a) 

where  Qn(x)  is  the  n ' th  Legendre  function  of  the  second 
kind,  it  can  be  shown  that 

A 

I(S)  =  2/  (2L+1)  { (L+l)  (L+2)  Ql+2(C)  -  L(L-l)  QL_2  ( C)  } 

-  (2/3)  <5L(1-  (IV- 3  lb) 


But 


Qn(x)  =  h  pn(x)  ln{  (x+1)  /  (x-1)  }  -  Wn_1(x), 

i-e.f  QR(x)  has  a  logarithmic  singularity  at  x  =  1.  in 
general,  numerical  integration  schemes  are  inaccurate  when 
such  behaviour  is  present,  and  one  has  to  extract  the 
singularity  analytically.  To  illustrate  the  principle, 
let  us  consider  integrals  of  the  type 


i 
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1  “  f  1  e  f(x>  1  g(y)/(x-y)  dy}  dx,  (IV-32a) 

where  f(x)  and  g(y)  are  well-behaved  functions,  for  which 
Taylor  expansions  about  unity  are  known: 


F(x)  =  ax  +  a2(x-l) 

+  a3(x-l)2 

-  f (x) , 

(IV-32b) 

G(y)  =  +  b2(l-y) 

+  b3(l-y)2 

-  g(y) . 

(IV-32c) 

Using  functions,  f(x)  and  g(y) ,  for  which  I  can  be 
evaluated  analytically,  extensive  calculations  showed  that 
the  following  procedure  works  extremely  well, 
i.  Evaluate 

I^x)  =  /_*  (g(y) -G(y)  }/(x-y)  dy  (IV-33a) 

numerically,  by  32— point  Gauss— Legendre  quadrature . ^ ® 
This  yields 

16 

I1(x)  *  Z  w  '  [(g(y.)-G(y.) }/ (x-y . )  + 

i=l  1  1  1  1 

(g(-yi) -G ( -y^ ) }/(x+yi) ] ,  (IV- 3 3b) 

where  w^  and  y^  are  the  Gauss— Legendre  quadrature 
weights  and  points.  Near  x  =  1,  I^x)  behaves  like 
(x-1)  In ( (x+1) / (x-1) } ,  and,  at  large  x,  like  1/x. 


i 


. 


(IV-34) 


ii.  Evaluate 

I2(x)  =  G(y)/(x-y)  dy 

=  b-j.  In {  (x+1) /  (x-1)  }  -  b2  [  (x-1)  ln{  (x+1)  /  (x-1)  }  -  2] 
+  b3  [(x-1)2  ln{  (x+1)  /  (x-1)  }  -  2  (x-1)  +  2] 

=  2  {blQo  -  b2(Ql-Q0)  +  (2b3/3)  (Q2-3Q1+2Q0)}. 

iii e  We  now  have 

1  ~  f(x)  (I1(x)+I2(x)  }  dx  (IV-35) 

=  ■/*!  e  x  (f(x)  I1(x)  +  [ f  (x)  -F  (x)  ]  1 2  ( x)  }  dx  + 

f°l  e"x  F(x)  I2(x)  dx. 

iv.  Evaluate 

I3  =  -ri  e  X  { f(x)  I1(x)  +  [  f  (x)  -F  (x)  ]  I2  (x)  }  dx  (IV-36) 

4  8 

by  Gauss-Laguerre  quadrature.  Near  x  =  1,  the 

term  in  braces  behaves  like  a  polynomial  plus 
3 

(x-1)  In { (x+1) / (x-1) } ,  and  the  numerical  integration 


is  very  accurate. 


* 


\  1  ,  .  -  . 


4 


Finally,  we  need 


I4  ^l  e  I2^x^  dx*  (IV-37) 

The  problem  here  is  the  logarithmic  singularity  at 
~  la  Direct  application  of  the  Gauss— Laguerre 
integration  formula  gives  very  low  accuracy.  How¬ 
ever,  it  is  possible  to  use  a  variation  of  this 
technique,  as  follows.  Consider 

N 

Tn  -  r°°  ~x  n  ,  ",  n 

I  -  /Q  e  x  dx  =  £  w  x.  ,  (IV-38) 

j  =  l  J  J 

where  x^  and  w^  are  the  quadrature  points  and 
weights.  If  n<N,  the  sum  is  exactly  In.  Now  set 

x  =  —  ln(t) ,  dt  =  -  e“x  dx. 

Equation  (IV-38)  then  becomes 
n  i  _  N 

1  =  ^*0  dt  =  E  w.  {-In  (t )  }n .  ( IV-39 ) 

j  =  l  J 

Again,  for  n<N,  the  sum  is  exactly  In ,  and  this 
gives  us  a  method  of  evaluating  I 


Write 


'•  • 

•• 
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where 

ti  =  1  +  exp(-xi)/  z±  =  2  +  Xi.  (IV-40b) 

We  now  return  to  the  problem  of  computing  the  matrix 
elements  of  (9Ve/3R)j  .  As  in  the  model  problem,  the 
basic  idea  is  to  isolate  terms  up  to  and  including  those 
which  behave  like  (5-1)  2  ln{  (  5+1)  /  ( 5-1)  }  ,  by  means  of 
Taylor  expansions  about  £  =  n  =  1. 

3-  Ia  =  <e  ’ r  ' I  OVe/3R)  *  | nFa>  (m=m'=0) 
a.  n-mtegrals 

Near  n  =  1,  expand  S^.^Cn)  S^(r))  as  a  power  series 
in  ( 1— rj )  .  Now 

PK°(n)  PL°(n)  ~  a0  +  ba(1_n)  +  ca(l-n)2,  (IV-41a) 


where 


bQ  =  “  h  (K  (K+l)  +L  (L+l)  }  ,  (IV-4 lb) 


Cq  =  d/16)  {  (K-1)K(K+1)  (K+2)  +  4K  (K+l)  L  (L+l)  + 


(L-l)  L(L+1)  (L+2)  }  .  (IV-4  lc) 


■ 


' 


Using  these  relations,  it  is  a  simple,  if  tedious,  matter 
to  write  (cf.  equation  (IV-29)) 
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=  XEQ(n)  +  xoQ(n) 


(IV-42) 


=  Aa  +  Ba(i-n)  +  cCT(i-n)2. 


From  equations  (IV-30)  and  (IV-31) ,  it  can  be  shown  that 


^anal 

O 


(Z)  -  U-o2)  -n  (?2-i)  }/U-n) 2  (a  +b  (i-n)  + 

GO 


ca(i-n)2}  dn 


(IV-43) 


=  4  {AaQ2  +  Ba(Q2-Q3)  +  (2Ca/15) (6Q4-15Q3+10Q2- 

Q0) }  +  (2/3) (Ba+2Ca) . 

Therefore,  Ia'(£)  in  equation  (IV-29)  is  given  by 

V(5)  =  Ud-n2)-n(C2-i)  )/(5-n)2  (xeo  (n)+xoQ  (n) }  dn 

=  i„nal  (5)  +  (C(l-n2)-n(C2-l)  }/(5-n)2  {xEo(n)  + 


xoa(n)-Aa-Bo(l-n)-CCT(l-n)2}  dn. 


(IV-44) 


-  * 


. 

< 
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The  leading  logarithmic  term  in  the  second  integral 
behaves  like  ( 5-1)  3ln{  (  5+1)  /  (5-1)  }  near  5=1,  and  I  used 
Gauss-Legendre  numerical  integration  to  evaluate  this 
integral  as  a  function  of  £;  call  the  result  The 

PL  s  can  ke  calculated  from  their  recursion  relation. 

b.  ^-integrals 

Equation  (IV-29c)  now  becomes 


I 

a 


(e2/4) 


'1 


{l^nal(5)+l"um(5) }  dC. 


Near  =1, 


(IV-45) 


xreg 

e  '  T  '  a 


(a  x 


nfa 


(C) 


exp(-cn5)  (5+l)q/2cn  1  {aa+6a  (C-l)  + 


Ya(C-l)2}, 


(IV-46a) 


where 


“a  “  a0b0 '  6o  =  alb0  +  h  aobl '  (IV-46b) 

Ya  =  *5  ( 2a2bo  +  aibi  +  ’5ao(b2~bl)  (IV-46c) 


and  a^ ,  b^  are  the  leading  terms  in  the  expansions,  (III- 
38)  (continuum  states)  ,  and  (III-17)  (bound  states)  ,  res¬ 
pectively.  Therefore, 


. 


■  ■ 


I 


a 


(e2/4)  /”  exp(-cnC)  (C)  F 


num 


6a(5-1)+YaU-l)2}  lanal(C)  de. 


(IV-47) 


As  in  the  model  problem,  the  first  integral  can  be 
performed  accurately  by  Gauss— Laguerre  quadrature;  and  the 
second  is  evaluated  in  two  parts:  from  1  to  2,  and  from  2 
to  oo  (See  equations  (IV-35)  to  (IV-40) ) .  The  Q  (£)  's  are 
computed  using  the  explicit  expressions  near  £= 1,  and  the 
power  series  at  larger  £ -values. 55 

4.  1^  =  <e,r,Tr|  (3Ve/3R)j  |nrTr>  (m'=m=l) 


g 


a„  ^-integrals 

We  follow  exactly  the  same  procedure  as  in  section  3 
above.  Near  n  =  1, 


PK1(n)  PL1(ri)  =  (1-rl)  ^a7r+b7I  ( 1-n ) } » 


( IV-4  8a) 


where 


=  %  K  (K+ 1)  L  (L+ 1) 


( IV-4  8b ) 


b 


7T 


h  a^  (K(K+1)  +  L  (L+ 1)  -  2}. 


(IV-4  8c) 


.  '  .  'N 
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Using  the  above  relations,  we  can  expand  the  product  of 
the  angular  eigenfunctions  as  a  power  series  in  (1-n) : 


q(sr  ,irsrTT+sr,wsriT)  +  r<s? 


=  XE^fn)  +  xo^tn) 


(IV-49) 


=  ( i— n ) 


{A  +B 

TT  TT 


(l-n) ) . 


From  equations  (IV-30)  and  (IV-31) , 

I*nal(5)  =  f+_l  (c (l-n2)  -n  (C2-i)  }/(C-n) 2  {^(i-n)  + 

2 

B^Cl-n)  }  dn  (iv-50 ) 

=  4  {A7t(Q2-Q3)  +  (2Bir/15)  (6Q4-15Q3+10Q2-Q0)}  + 

(2/3) (A  +2B  ) ; 

TT  7T 


and 


.num 


.+  1 


I7r  (?)  =  f_1  Ud-n  )-n(€  -i)}/(C-n)  (xe^  (n)  txo^  (n)  - 


( l-n) -b^  ( l-n)  }  dn 


+  1  16 

H(n)  dn  -  £  w.  '  {h (n  .  )+H(-n . ) } , 

i=l  1  1  1 


i 


(IV-51) 


r 


•• 


where  r^,  w^'  are  the  Gauss -Legendre  quadrature  points  and 
weights . 
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b.  ^-integrals 

I7r  =  (e2/4)  /“  exp(-cnC) 


xreg 

e  '  r 


7T 


(5) 


W«> 


{lanal(5) 


j.num 

7T 


a  >  dc. 


( IV-52 ) 


Near  £  =  1, 


xreg 
e  T 


TT 


(O 


Fnr,(«) 


=  (5-1)  (5+l)q/2cn_2{air+B7T(C-l)  }, 


where 


(IV-53a) 


aTT  “  a0b0  '  3tt  =  alb0  +  *  a0bl7  ( IV-53b) 

and  a^ ,  are  the  leading  terms  in  the  expansions,  (HI- 
38)  (continuum  states)  ,  and  (III-17)  (bound  states)  ,  res¬ 
pectively.  Therefore, 

Xir  =  (e2/4)  fl  exP(-cn5){Xe'r'^(5)  Fnr7r(?)  + 

(xe*r-ir(5)  Fnrw(«>  "  (5-D  (5+1) q/2cn"2  [c^+B, «-l)  ]) 
Ianal(5))  <35  +  (e2/4)  /”  exp (-c  5)  ( 5-1)  ( 5+1)  q,/2cn~2  {a  + 

x  n  77 

6^(5-!)}  Ianal  ( 5)  d?. 


(IV-54) 


. 


. 


•• 


•  • 


110 


These  integrals  are  computed  exactly  as  in  the  I 


case . 


This  completes  the  calculation  of  the  matrix  elements 

given 


of  OVe/9R)j  .  The  radial  component  of  P(R)  is  then 
by  9 


V<V£i>- 


(IV-55) 


We  now  turn  to  the  matrix  elements  of  sg. 


E.  Matrix  Elements  of  \h _ sgl 

■  — - BO - 


1.  Introduction 

We  have  to  compute  (see  equations  (IV-5b)  and  (IV- 


9b)) 


I  R  =  <e,r,m|z  f | nTm> 

m  g  1 


(IV-56a) 


zg  =  (R/2) 


(IV-56b) 


ii-  Im®m•  =  |xgf |nrm> 


(IV-57a) 


xg  =  (R/2)  (C2-l)!5(l-n2)*5 


COS(J) 


( IV-5  7b) 


For  the  switching  function,  I  propose  to  use  the 
form  (see  figure  7) : 


f  =  tanh ( An+B) , 


( IV-5  8) 


•  ^ 


. 


. 


Ill 


Figure  7.  Model  switching  function:  f  =  tanh ( 2 n  —  0 . 2) 


. 
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where  the  parameters,  A  and  B,  are  to  be  determined.  This 
will  be  discussed  in  detail  in  the  next  chapter. 

The  integrals  over  (p  can  be  done  by  inspection,  and 
yield  the  selection  rules  for  m:  matrix  elements  are  non¬ 
zero  only  if  m1  =  m,  for  the  radial  component;  and  m'  =  m 
-  1,  for  the  angular  component.  Of  course,  these  rules 

are  the  same  as  those  for  the  matrix  elements  of  ( 9V  /9R)->- 
~  e'  r 

and  L±g,  respectively.  Both  the  n~  and  ^-integrals  are  ^ 

evaluated  numerically,  and  they  are  quite  straightforward. 

2.  Evaluation  of  I  R 

m 


!mR  =  (R/2)4  / 


1  e*P<-cnS>  Xe'?'i „<«>  W«  5  dC, 

(IV-59a) 


where 


xl  "  ^-1  sr'm(n)  Srm(n)  n  tanh(An+B)  dn,  (IV-59b) 
+  1  3 

I2  =  •'-1  sr'm^  Srm^n)  n  tanh  (Ar>+B)  dn.  ( IV-59  c) 


The  integrals,  1^  and  ,  are  evaluated  by  32-point 
C^uss -Legendre  quadrature.  To  simplify  the  expressions, 
write  n  , m ( n )  Spm(n)  as  the  sum  of  even  and  odd  parts. 


n  srim(n)  srm(n)  =  SE(n)  +  so(n). 


(IV-60) 


which  can  be  calculated  explicitly  from  equations  (IV-13a) 


‘  «q  yi  1  M  1 


<  ,  &U  tr** 


. 


♦ 
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and  (IV-15a) .  Therefore, 


16 


1  ~  Z  w-i  '  {[SE(n.)+SO(n.)]  tanh  (An  .  +B)  -  [SE  (n  . ) -SO  (n  . )  ] 
j=l  J  J  J  j  33 


tanh (-An j+B) };  (IV-61a) 


and 


16 


I2  ~  1  wj  '  t  [SE(r]j)+SO(ri  j)  ]  tanh(Arij+B)  -  [SE  (n  j  )  -SO  (n_.  )  ] 


tanh (-An j+B) }  r |_.  .  (IV-61b) 


The  ^-integrals  are  now  evaluated  by  Gauss-Laguerre 
quadrature . 


3.  Evaluation  of  I 


0 

m+m 


4  -OO 


Vm'  =  15  (R/2)’  exp(-cnS)  Fnrm(5)  (S2-l) 


2 


U  I3-I4)  d$. 


(IV-62a) 


where 


X3  -  f+i  ST,m, (n)  srm(n)  (l-n2)^  tanh(An+B)  dn ; 


( IV-62b) 


and 


T 


■ 


. 


I4  ~  ^-l  Sr'm'^^  Srm^r'^  (1_n2)  ^  tanh(Ar|+B)  r]2  dn . 

(IV-62c) 

The  factor  of  %  comes  from  the  integration  over  <p.  Now^^ 

U-n2)*5  PLm(n)  =  (2l+i) -1  {PL^1(n)-PL”;1(n)}, 

and  this  relation  can  be  used  to  simplify  the  integrals, 
so  that  they  only  involve  one  set  of  P  m ( n )  's.  If  m'  =  m 

Li 

+  lf  use  the  relation  to  express  (l-p2)^  Sr  (n)  as  a  sum 

rm 

in 1 

ef  ( ^ )  s;  on  the  other  hand,  if  m'  =  m  -  1,  use  it  to 

express  (l-n2)^  sr.m.(n)  as  a  sum  of  PLm(n)'s.  In  either 
case,  write  ( 1— rj  )  ^  Sp,m,(n)  S^Cri)  as  a  sum  of  even  and 
odd  parts : 

(i-n2P  sr,m,(n)  srm(n)  =  ye (n)  +  Yo(n).  (iv-63) 

Thus,  I ^  and  1^  reduce  to  exactly  the  same  forms  as  1^ 
(equation  (IV-59b) )  and  I2  ((IV-59c)),  and  the  evaluation 
of  proceeds  on  exactly  the  same  lines  as  that  of  I  R. 

The  components  of  A ( R)  are  then  given  by 

tA9,jn  =  <im/2ti)  (ej_en)  ImR'  (IV-64a) 

[Ag]®n  =  (im/2h)  (£.-en)  (IV-64b) 


'■ 


. 


-  5  f  ;  *  : 


•  ■. 


' 


115 


F.  Summary 

In  this  chapter,  I  have  described  the  evaluation  of 

the  matrix  elements  af  P(R)  +  A(R)  in  great  detail,  for 

transitions  from  bound  molecular  states  to  the  electronic 

continuum.  For  the  next  chapter,  I  would  like  the  reader 

to  imagine  that  we  have  a  "black-box"  computer  programme. 

If  we  give  the  programme: 

i.  the  nuclear  charges,  e.g.,  Z  =  2.0,  Z  =  1.0 

A  B 

ii.  bound  state,  united-atom-limit  quantum  numbers,  e.g., 
n=2,  F=l,  m=0  (2pa  state) 

iii.  continuum  state  energy  (or  energies) ,  united-atom- 
limit  angular  momentum,  and  azimuthal  quantum 
numbers,  e.g.,  e ’  =  0.5  Ry. ,  T'  =  1(1)5,  m'  =  1 
(p,  d,  f,  g,  and  h  7r-waves) 

iv.  the  parameters,  A  and  B,  in  the  switching  function 

v.  the  desired  range  of  R-values,  e.g.,  0.25(0.25)10.0 
then  the  computer  will  calculate  for  us,  at  each  R-value: 
i.  the  bound  state  eigenenergy  and  separation  constant 
ii.  the  continuum  state  separation  constant  and  phase 
shift 

iii.  [Pg(R)]  and  [Pg (R) +Ag (R) ] . 

Dn  jn 

In  our  example,  it  will  calculate  the  angular  couplings 
from  the  2pa  state  to  the  five  partial  waves  specified. 

In  the  next  chapter,  I  will  discuss  the  variation  of  the 
matrix  elements  with  the  parameters,  A  and  B. 
assume  that  the  programme  is  debugged. . . 


We  will 


. 


»  • 
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CHAPTER  V 


RESULTS 

A.  Introduction ^ 

In  earlier  chapters ,  I  have  pointed  out  that  non— 
adiabatic  couplings  between  low-lying  bound  states  and  the 
electronic  continuum  are  extremely  sensitive  to  the  choice 
of  the  switching  function.  Studies  of  this  f-dependence 
may  provide  useful  insight;  and  it  may  also  be  possible  to 
choose  an  f  which  is  in  some  sense  "optimum"  for  a  partic¬ 
ular  problem. 

Such  studies  have  been  carried  out  by  Thorson  and 

coworkers,2  in  connection  with  the  calculation  of  direct 

impact  ionisation  cross  sections  in  H+-H(ls)  collisions. 

For  each  bound  state  they  studied  (Isa  ,  2pa  ,  2dtt  ,  in 

g  u  u 

united  atom  notation) ,  they  found  that  a  definite  choice 
for  the  switching  function  may  be  made,  which  simulta¬ 
neously  and  systematically  reduces  the  couplings  from  that 
state  to  all  continuum  states,  in  most  cases  by  several 
orders  of  magnitude.  Whereas  the  uncorrected  pss  theory 
predicts  large  couplings  to  30-40  continuum  partial  waves, 
with  the  envelope  of  these  couplings  having  a  range  of 
about  40  au;  the  corrected  couplings  are  significant  only 
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for  the  first  two  or  three  partial  waves,  and  their  range 
is  less  than  12  au.  It  was  shown  that  the  switching 
functions  which  achieve  these  spectacular  reductions  do 
not  depend  on  the  energy,  or  other  quantum  numbers,  of  the 
continuum  states,  nor  upon  the  type  of  coupling  causing 
transitions  (radial  or  angular);  however,  the  f's  are 
different  for  each  discrete  state.  From  these  results, 
one  could  infer  that  the  switching  functions  so  determined 
are  characteristic  of  the  discrete  states  themselves,  and 

and  might  therefore  be  relevant  to  the  computation  of 
corrected  couplings  for  transitions  other  than  ionisation 
for  example,  close  coupling  between  discrete  states. 

This  possibility  was  an  important  motivation  for  the  work 
reported  in  this  chapter.  I  have  extended  the  studies 
described  in  reference  2  to  electrically  asymmetric  sys¬ 
tems,  such  as  HeH  ,  and  to  additional  discrete  states;  I 
have  also  made  the  procedure  for  selecting  the  switching 
function  parameters  more  systematic.  The  corrected  matrix 
elements  reported  here  can  also  be  used  to  calculate  ioni¬ 
sation  cross  sections,  for  these  more  general  systems. 

The  formulation  of  slow  collision  theory  using 
switching  functions  can  be  extended  so  that  each  discrete 
electronic  state,  |n>,  has  a  translation  factor  with  a 
different  characteristic  switchina  function,  f 

n 

This  extension  is  quite  straightforward  at  the  level  of 
semiclassical  approximation,  and  the  main  results  are  as 


a 
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follows  (see  Appendix) . 

i.  The  resulting  close-coupled  equations  have  the  same 
form  as  those  derived  in  chapter  II ,  and  the 
expression  for  the  corrected  nonadiabatic  coupling 

ma*"r:*-x  elements  also  has  the  same  form  (see  equations 
(V-2)  and  (V-3) ,  below) 

11.  Since  the  switching  functions,  f n ,  f  ,  are  different, 
the  electronic  states,  [n>,  [j>,  are  not  orthogonal, 
the  corrected  nonadiabatic  couplings  are  not 
Hermitian;  however,  the  non— Hermitian  character 
exactly  compensates  for  the  non-orthogonality  of  the 
basis  kets ,  and  guarantees  that  the  coupled  equations 
conserve  probability 

In  the  asymptotic  channel  limits,  the  electronic 
states  are  rigorously  orthogonal,  so  vector  solutions 
to  the  coupled  equations  give  directly  the  unitary 
S-matrix  for  a  collision,  in  the  limit,  t  +°° 

This  establishes  that  a  switching  function  formu¬ 
lation  of  a  close-coupling  slow  collision  problem,  with  a 
different  switching  function  for  each  molecular  state, 
provides  a  formally  valid  theory  (provided,  of  course, 
that  the  f  ' s  satisfy  the  asymptotic  constraints). 

As  I  pointed  out  in  chapter  I,  most  earlier  treat¬ 
ments  of  the  slow  collision  problem,  which  correctly  take 
account  of  boundary  conditions  and  translation  factors, 
can  be  recovered  as  special  cases  of  this  general 


' 
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formulation.  In  electrically  asymmetric  systems,  the 

molecular  states  asymptotically  correlate  in  a  one-to-one 

fashion  with  atomic  states  of  either  nucleus  A  or  B.  If 

we  choose  f  =  +1  for  the  asymptotic  B-states,  and  f  =  -1 

n 

for  the  A-states,  the  formulation  of  Bates  and  McCarroll3 

is  obtained.  This  is  also  exactly  equivalent  to  the 

formalism  (based  on  projection  operators)  proposed  by 
2  ^ 

Matveyenko,  and  is  essentially  equavalent  to  the  "dis¬ 
torted  cluster  states"  method  of  Y.  Hahn.23  For  symmetric 
systems,  the  connection  between  a  switching  function 
approach  and  these  other  methods  is  not  so  clear. 

In  this  chapter,  I  will  describe  the  selection  of  a 
more  elaborate  set  of  switching  functions  for  several 
discrete  states  of  two-centre,  one-electron  systems.  The 
states  studied  are  the  Isa,  2sa ,  2pa ,  2pTr,  and  3da  states 
(united  atom  notation);  for  systems  with  nuclear  charges, 

~  +1.0,  +1.0  <  ZA  <  +2.0.  The  selection  criterion  is 
based  on  the  remarkable  sensitivity,  of  continuum  coup¬ 
lings  from  these  states,  to  the  detailed  form  of  the 
switching  function,  as  reported  earlier  for  H  +.  We  find 
that  there  are  very  specific  choices  for  each  switching 
function,  f^(r;R) ,  which  are  independent  of  continuum 
state  properties,  or  the  type  of  coupling  considered,  but 
are  characteristic  of  each  discrete  state,  |n>.  As  stated 
earlier,  this  approach  suggests  that  these  switching 
functions  may  be  appropriate  for  problems  more  general 


. 
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than  the  ionisation  one. 


These  calculations  are  not  based  rigorously  upon  any 
formal  variational  principle,  nor  is  the  set  of  switching 


functions  so  obtained  unique,  or  a  "best,"  choice  in  any 
absolute  sense.  However,  this  set  may  have  some  advan¬ 
tages  for  close-coupling  studies,  and  this  point  is  dis¬ 
cussed  in  section  D,  by  some  comparisons  with  alternative 
formulations . 

In  the  rest  of  the  Introduction,  I  will  summarise 
the  definitions  of  the  coupling  matrix  elements.  Section 
B  describes  the  procedure  used  to  select  the  switching 
function  parameters;  and  section  C  presents  the  results  of 
the  calculations. 

In  un corrected  pss  theory,  the  matrix  elements  for 
nonadiabatic  coupling  are  defined 


(V-l) 


In  the  corrected  theory,  P(R)  is  simply  replaced  by  Pg(R) 
+  Ag ( R) ,  where 


( V-2 ) 


and 


( V-3) 


' 
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with 


s 

n 


9  =  fn(?;R) 


(V-4) 


These  equations  are  exactly  equivalent  to  equations  (IV-1) 
to  (IV-4) ,  with  f (r ; R)  replaced  by  fn(r;R)  (see  Appendix). 

B.  "Optimisation "  Procedure 

For  their  calculations  on  the  H2+  system,  Lebeda, 

2 

Thorson,  and  Levy  a  used  a  switching  function  of  the  form 


fn(r;R)  =  tanh{Bn(R)Rn>,  (V-5) 

where  3r(R)  is  a  variable  parameter.  For  each  discrete 
state  they  studied,  at  a  series  of  R-values,  they  plotted 
the  magnitudes  of  the  corrected  radial  coupling  matrix 
elements  against  3,  for  a  series  of  partial  waves.  At 
each  R,  they  showed  that  there  is  a  particular  value,  3 

n ' 

which  produces  spectacular  reductions  in  the  corrected 
couplings  to  higher  partial  waves,  relative  to  the  pss 
values.  The  same  value  of  3r  was  obtained  by  independent 
calculations  at  several  different  continuum  energies. 
Similar  calculations  by  SethuRaman,  Thorson,  and  Lebeda2b 
for  the  angular  couplings  gave  approximately  the  same 
values  of  Br  as  were  obtained  for  radial  couplings. 

I  have  repeated  these  calculations,  and  extended 
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them  to  additional  discrete  states,  and  to  electrically 
asymmetric  systems;  rather  than  determine  Bn(R)  graphic¬ 
ally,  I  have  made  the  selection  procedure  automatic,  and 
numerically  more  precise.  However,  the  main  idea,  and  the 
behaviour  observed,  are  the  same  as  in  reference  2.  Using 
a  suitable  parametric  form  for  the  switching  function,  we 
look  for  specific  choices  of  the  defining  parameters  which 
achieve  large,  simultaneous  reductions  in  as  many  of  the 
continuum  couplings  as  possible,  and  we  define  as  "optimum" 
the  parameters  which  best  achieve  this.  This  idea  is 
justified  by  the  remarkable  results  obtained,  rather  than 
by  any  a  priori  arguments. 

For  a  system  with  nuclear  charges,  Z  ,  Z  ,  at  each 

H  D 

R,  and  for  each  discrete  state,  Jn>,  considered,  the  pro¬ 
cedure  that  I  used  is  as  follows. 

i.  Consider  a  set  of  continuum  states,  { <  j  |  } ,  with  the 
same  (arbitrary)  energy,  e  ^  ,  and  azimuthal  quantum 
number,  m';  and  linked  to  |n>  by  either  radial 
(m'=m)  ,  or  angular  (m'=m±l)  ,  coupling.  The  set, 

( < j I  ) ,  is  thus  a  set  of  "partial  waves"  characterised 
by  the  "angular"  eigenvalues,  and  I  have  indexed 
them  by  the  corresponding  united-atom-limit  orbital 
angular  momentum  quantum  number,  L' . 
ii.  For  each  state,  <j|,  we  define  the  residuum 
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where  the  matrix  elements  are  the  radial  (equation 
(IV-5) ) ,  or  angular  (equation  (IV-9)),  components  of 
—  and  A  (with  f  replaced  by  f^)  .  Pjn  the 

magnitude  of  the  corrected  coupling,  relative  to  the 
uncorrected  (geometric  centre)  pss  matrix  element. 

iii.  The  parametric  form  used  for  f  is 

n 

fn  =  tanh{R(!5Bn[(ZA+ZB)n+(ZA  ZB)  ]  +  anln  (ZB/ZA)  )  }  (V-7) 


where  an  and  $n  are  variable  parameters.  This 
reduces  to  the  form  ( V  5 ) ,  for  Z^  =  Z^;  it  is  equiv¬ 
alent  to  the  form,  tanh  {R  (B^ri+A^)  }  (cf.  equation  (V- 
58) ) ,  but  it  displays  more  explicitly  the  observed 
charge  dependence  in  the  optimum  switching  functions . 
iv.  For  each  set  of  continuum  states,  we  determine  the 
parameters,  an  and  3  ,  such  that  the  sum 


S  (a  ,6  )  =  Z  p. 
n  n  n  .  in 

J  J 


(V-8) 


is  minimised.  Most  of  the  ratios,  PjR,  are  very 
sensitive  to  changes  in  f  :  the  residua  of  the  first 
two  or  three  partial  wave  couplings  are  only  changed 
slightly  by  varying  f  ;  but  those  for  all  higher 
partial  waves  vary  by  orders  of  magnitude,  and  the 
same  choice  of  f  reduces  all  of  them  together.  The 
pjn  are  tYPically  much  more  sensitive  to  3  than  to  a. 
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The  sum  in  equation  (V-8)  may  include  all  of  the 
partial  waves  considered  (up  to  L*  =  10  or  12),  or  it  may 
include  only  some  of  them— as  long  as  the  higher  ones  are 
included.  In  searching  for  the  minimum,  I  found  it  con¬ 
venient  to  omit  the  first  two  or  three  partial  wave  coup¬ 
lings,  because  of  their  insensitivity  to  changes  in  f  ; 
the  minimum  is  much  more  sharply  defined  when  they  are 
excluded.  However,  the  "optimum"  values  obtained  for 
8n(R)  (to  three  significant  figures),  and  an(R)  (to  two 
sig.  figs.)  remain  unchanged. 

This  process  is  repeated  for  other,  arbitrarily 
chosen,  energies,  and  for  all  possible  types  of  coupling 
from  the  state,  |n>,  at  a  series  of  R-values;  and  the 
values  of  the  "optimum"  parameters  are  compared.  The 
results  are  described  in  the  next  section. 

C.  Results  of  Calculations 

1.  Behaviour  of  Optimum  Parameters 

(a)  The  precision  with  which  the  optimum  parameters 
can  be  determined  depends  on  three  factors:  (1)  the  sharp¬ 
ness  of  the  minima  in  S^,  (2)  the  agreement  between  deter¬ 

minations  using  different  continuum  energies,  and  (3)  the 
agreement  using  different  types  of  coupling.  Usually, 
these  factors  are  commensurable,  and  we  can  determine  6 

n 

to  three  significant  figures  and  to  two  significant 
figures . 
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(b)  In  every  case,  the  optimum  parameters  are 
independent  of  continuum  energy. 

(c)  In  all  cases  (except  for  the  3da  state;  see 
items  (e)  and  (h) ,  below) ,  the  optimum  parameters  are  the 
same  for  radial  couplings  from  |n>  as  for  angular  coup¬ 
lings;  in  particular,  this  is  true  for  the  2pr  state, 
which  is  linked  by  angular  couplings  to  both  a  and  5 
continuum  states. 

(d)  I  have  redetermined  the  parameters,  3  (r)  for 

n 

both  radial  and  angular  couplings  from  the  Isa  ,  2pa  ,  and 

|  ^  ^ 

2p1Tu  states  of  h2  '  which  were  reported  in  reference  2. 

The  new  values  lie  well  within  the  (graphically  estimated) 
ranges  given  in  reference  2a,  for  the  radial  parameters, 
but  do  not  always  agree  with  the  estimated  angular  param¬ 
eters  of  reference  2b.  The  present  numbers  are  the  more 
reliable. 

(e)  I  have  also  determined  3  (R)  for  the  2sa  and 

+  n  g 
3dCg  states  of  H 2  .  However,  for  the  3da  state,  the 

9 

angular  couplings  are  not  very  sensitive  to  the  choice  of 
3r  and  the  resulting  minimum  is  nearly  flat  for  0.02  <  3  < 

0.10.  I  have  taken  3R0  =  3rR  here,  as  is  found  for  every 
other  case. 

(f)  I  have  studied  the  Isa,  2s a,  2pa,  2p tt,  and  3da 

states  of  the  asymmetric  HeH++  system,  and  the  Isa,  2p a, 

and  2 pTr  states  of  the  intermediate  model  systems  with  Z  = 

B 

1.0,  and  —  1.2,  1.4,  1.6,  and  1.8.  If  f  is  described 
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by  equation  (V-7) ,  the  effects  of  charge  asymmetry  are 
very  well  accounted  for,  and  we  can  take  the  optimum  values 
°f  &n(R)  to  be  the  same  as  those  for  H2+.  Typically,  the 
Sn  surface  has  a  deep,  somewhat  elongated  depression  in 
the  (a, 3)  plane,  running  almost  parallel  to  the  a-axis. 
Within  this  entire  region,  the  residua  for  the  higher 
partial  waves  are  all  of  the  order  10~3  to  10~4.  There- 
^ore '  — rc}al b  changes  in  3^(R)  may  be  offset  by  compensating 
j££9er  changes  in  an(R),  while  remaining  in  the  optimum 
region.  In  most  cases  (but  see  item  (g) ,  below),  the 
range  of  acceptable  3-values  does  not  vary  from  the  H  + 
values  by  more  than  1%  to  4%.  The  important  point  is  that 
no  significant  further  reduction  of  coupling  matrix 
elements  occurs  elsewhere  in  the  valley. 

The  optimum  ( R)  vary  slowly  with  charge  asymmetry. 

(g)  At  distances,  R  >  2  to  4  au,  the  optimum  param- 
eters  for  the  lscr  and  2sg  states  of  HeH"*""*"  cannot  be  so 
precisely  determined.  Even  at  short  distances,  these 
states  are  dominantly  atomic  states  of  the  He+  ion,  and 
the  appropriate  switching  function  is  -1.  There  is  a  wide 
range  of  parameters,  a  and  3,  yielding  that  simple  result. 

larger  R— values  (R  ^  12  au) ,  this  becomes  the  charac¬ 
teristic  situation  for  most  of  the  molecular  states  of 
HeH++;  in  effect,  the  deep,  elongated  valley  described 
above  becomes  very  long  and  very  broad.  This  is  to  be 
expected,  because  the  very  simple  "atomic"  switching 


. 


4 


-  '  ... 


.  '  - 


4 


-  } 

• 

V 

* 

» 

...  . 

- 

. 


functions ,  f  —  fj,  are  always  asymptotically  correct. 
However,  the  Isa  and  2sa  states  of  HeH++  are  the  only 
cases  where  such  very  simple  choices  can  produce  efficient 
reduction  of  the  spurious  pss  couplings  at  the  small  to 

moderate  distances  where  the  physically  real  couplings 
occur. 

(h)  As  in  the  case  of  H2  ,  I  found  that  the  angular 
coupling  matrix  elements  for  the  3da  state  of  HeH++  are 
not  very  sensitive  to  changes  in  a  and  8,  and  I  have  again 
taken  a0  =  aR,  and  3°  =  3R. 

Figure  8  shows  the  values  of  8n(R)  vs.  R,  determined 
for  .  For  the  3da^  state,  the  crosses  represent  the 
values  determined  for  the  radial  couplings  alone.  For  the 
remaining  states,  where  an  error  bar  appears,  it  indicates 
the  difference  between  "best"  values  for  3  ( R)  obtained 
from  radial  and  from  angular  couplings;  where  only  a  cross 
appears,  this  difference  is  too  small  to  be  shown.  Figure 
9  depicts  a^CR)  vs.  R,  for  HeH++  (8  values  are  those  for 
h2+) .  Error  bars  and  crosses  have  the  same  meaning  as  in 
figure  8.  The  large  uncertainties  in  alga  and  ct2sa  are 
related  to  the  simple  behaviour  of  the  switching  functions 
for  these  states. 

Figures  10  and  11  are  plots  of  the  "optimum" 
switching  functions  for  H2+  and  HeH++  vs.  n,  at  R  =  8.0  au. 
In  figure  12 ,  I  have  shown  the  behaviour  of  the  switching 
function  for  the  Isa  state,  as  a  function  of  Z  ,  again  at 
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Figure  8.  "Optimum"  0^'s,  for 


(au) 
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Figure  9.  "Optimum"  an*s ,  for  HeII++ 


(au) 
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Figure  10.  "Optimum"  f  's,  for  H  + 

n  2  ' 


at  R 


8.0  au 
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Figure  11.  "Optimum"  f^'s. 


++ 


for  HeH 


,  at  R  =  8.0  au 
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Figure  12.  "Optimum"  ^  *s,  as  a  function  of  Z  , 

A 

with  Z0  =  1.0,  and  R  =  8.0  au 
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R  =  8.0  au.  This  illustrates  very  clearly  the  development 
of  the  atomic  character  of  f^,  as  ZA  increases. 

Table  5  gives  detailed  results  for  the  "optimum" 
parameter,  6n#  Tables  6  to  9  show  the  "optimum"  values  of 

cl  ,  as  a  function  of  Z,  . 
n  A 

2.  Matrix  Elements 

To  compute  corrected  nonadiabatic  couplings  to  the 
continuum,  I  have  used  the  straight-line  approximations 

k°th  and  an  (R)  vs.  R,  which  are  shown  in  figures 

8  and  9.  It  would  appear  from  the  figures  that  this  might 
introduce  large  errors  in  the  matrix  elements  at  small 
values  of  R,  but  this  is  not  the  case.  For  R  <  4  au,  the 
matrix  elements  to  higher  partial  waves  are  extremely 
small,  while  the  couplings  to  the  lower  partial  waves  are 
not  very  sensitive  to  changes  in  the  switching  function 
parameters.  Whether  such  departures  from  the  "optimum" 
parameters  can  be  used  in  computing  discrete  to  discrete 
couplings  is  not  so  clear. 

I  have  calculated  the  corrected  matrix  elements  for 
the  five  discrete  states  mentioned,  for  H2+  and  HeH++. 

The  corrected  matrix  elements  for  all  but  the  first  few 
partial  waves  are  negligible,  and  these  few  are  usually 
significantly  reduced  in  size  and  range.  By  contrast,  the 
pss  couplings  are  all  large  and  very  long  ranged. 

Notation  in  the  figures  labels  bound  states  (|nLm>), 
and  continuum  states  ( <e ' L 1 m ' | ) ,  where  n  is  the  united 
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fl 

OPTIMUM" 

TABLE  5 

3n's,  DETERMINED  FOR 

H2+ 

Isa 

2sa 

2pa 

R 

A 

L 

y 

9/9R 

Ly  3/ 3K 

A 

L 

y 

9/9R 

2 

.  399 

.400 

.2231  .2239 

.212 

.203 

4 

.  385 

.  385 

.2050  .2053 

.216 

.217 

6 

.402 

.403 

.2033  .2037 

.220 

.220 

8 

.424 

.424 

.2073  .2073 

.224 

.223 

10 

.439 

.439 

.2136  .2135 

.229 

.227 

12 

.450 

.451 

.2200  .2201 

.237 

.231 

2pir 

3da 

A 

A 

R 

L_ 

L+  9/9R 

9/  9R 

2 

— 

.161 

.042 

4 

e  161 

.160  .152 

.063 

6 

.  157 

.155  .146 

.063 

8 

.  149 

.151  .151 

.058 

10 

.  152 

.154  .154 

.052 

12 

.  159 

.159  .159 

.046 
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TABLE  9 


"OPTIMUM" 

Vs' 

FOR  HeH**"1' 

(a) 

2sa  state 

/s 

R  =  2 

4 

6 

8 

10 

12 

L 

y 

0.88 

1.43 

1.05 

0.95 

1.2 

1.2 

3/3R 

1.6 

1.28 

0.91 

0.92 

1.0 

0.9 

(b)  3d a  state 

R  -  2 

4 

6 

8 

10 

12 

3/3R 

0.55 

0.23 

0.17 

0.17 

0.19 

0.24 

•• 
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atom  principal  quantum  number,  e1  is  the  continuum  elec¬ 
tronic  energy  (in  Rydbergs);  L,  L1  are  united  atom  orbital 
angular  momentum  quantum  numbers;  and  m,  m'  are  azimuthal 
quantum  numbers;  distance,  R,  is  in  au.  The  matrix 
elements  shown  are  related  to  those  defined  in  equations 
(IV— 5)  and  (IV— 9)  as  follows: 

<e 'L 'mf  H  (RAD)  |  nLm>  =  (i/fi)  [Pg+Ag]*n,  (V-9a) 

<£ 'L ' m± 1 [ H (ANG) | nLm>  =  (i/2yfi)  [Pg+Ag]°n.  (v-9b) 

For  convenience,  I  have  taken  y  =  1836.1  mQ  (proton  rest 
mass).  Table  10  provides  a  guide  to  the  figures. 

Figure  13  shows  the  disastrous  long-range  behaviour 
of  the  uncorrected  pss  matrix  elements,  here  shown  for 
H2  7  the  corresponding  corrected  matrix  elements  (dashed 
curves)  are  negligible.  Of  course,  at  the  large  distances 
shown,  any  simple  switching  function  which  has  the  value 
+1,  in  the  entire  neighbourhood  of  nucleus  B,  and  -1,  near 
A,  will  produce  similar  results;  the  main  point  of  this 
picture  is  to  illustrate  the  uncorrected  pss  couplings. 

Figures  14  to  24  compare  the  uncorrected  pss  coup¬ 
lings  (solid  curves)  and  the  corrected  couplings  (dashed 
curves),  for  the  various  initial  states  of  H  +;  figures  25 
to  35  show  the  same  thin  or  HeH++.  For  HeH++,  I  have 
also  included  for  compare,  .m  the  corrected  matrix  elements 
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TABLE  10 

A  GUIDE  TO  THE  FIGURES 


Isa 

2sa 

2pa 

2p7T 

3da 

H2  + 

14 

15 

16 

17 

18 

radial 

HeH++ 

25 

26 

27 

28 

29 

H2  + 

19 

20 

21 

22  (LJ 

24 

angular 

23(L+) 

HeH++ 

30 

31 

32 

33 (L_) 

35 

34 (L+) 

- - —  pss  value 

corrected  value  with  optimum  f 

n 

corrected  value  with  Bates-McCarroll  f 

n 

In  all  cases,  the  continuum  energy,  e',  is  0.50  Ry. 
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Figure  13.  Long  range  behaviour  of  pss  matrix  elements, 

shown  here  for  H  * 
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Figure  14.  <0 . 5L ' 0 | H ( PAD)  j 10 0 > ,  H2+,  L'  =  0(2)8 
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Figure  16. 


<0. 5L' 0 |H (RAD) | 210> ,  H2+,  L' 


1(2)  7 
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Figure  17.  <0 . 5L ' 1 | H (RAD) | 2 11> ,  H„+,  L'  =  1(2)7 
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Figure  18.  <0 . 5L '  0  |  H  (RAD)  |  320>  ,  H2+,  L' 


0(2)8 


. 


147 


.. 


•s  -  ' 


148 


8 


CO 


CO 


C\J 


Figure  20. 


<0.5L'l|H(ANG)  |200>,  H  +,  L' 


=  2(2)8 
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Figure  21. 
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Figure  22.  <0 . 5L  '  0  |  H  (ANG)  |  2 11>  ,  H2+,  L'  =  1(2)7 
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Figure  23. 
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Figure  24.  <0 . 5L ' 1 | H (ANG) | 320> ,  H2+ ,  L'  =  2(2)8 
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Figure  28.  <0 . 5L ' 1 | H (RAD) | 2 11> ,  HeH++,  L' 
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Figure  30.  <0 . 5L ' 1 | H ( ANG)  j 100> ,  HeH 
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Figure  31. 
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Figure  32. 
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Figure  33.  <0 . 5L 1 0 | H ( ANG)  | 2 11> ,  HeH++ ,  L'  =  0(1)5 
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Figure  34 . 
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Figure  35.  <0 . 5L ' 1 1 H ( ANG)  | 320>  ,  HeH++ ,  L'  =  1(1)6 
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which  are  obtained  using  the  simple,  Bates-McCarroll 3 
switching  functions  (fR  =  +1,  for  asymptotic  "B"  states; 

^n  ~~  ^or  ^  states)  .  As  I  noted  above,  this 

simple  choice  gives  excellent  results  for  the  Isa  and  2sa 
states  (figures  25  and  26)  at  all  distances.  For  the 
other  states,  this  choice  gives  results  in  agreement  with 
those  using  "optimum"  switching  functions,  at  large 
distances  only.  For  intermediate  distances,  the  resulting 
couplings  are  often  larger  than  the  uncorrected  pss  coup¬ 
lings.  This  is  especially  true  for  states  like  2pa,  which 
have  significant  molecular  character  (figure  27).  These 
comparisons  show  that,  in  such  cases,  the  more  elaborate 
description  of  switching  functions  given  here  is  really 
necessary  to  describe  continuum  couplings  well.  This  will 

probably  also  be  true  for  couplings  between  discrete 
states . 


D.  Discussion 

1.  Summary 

Using  switching  functions  to  describe  the  effects  of 
electron  translation  in  discrete  molecular  states,  I  have 
shown  empirically  that  (1)  the  corrected  nonadiabatic 
couplings  from  these  states  to  the  continuum  are  generally 
very  sensitive  to  the  detailed  form  of  the  switching 
function;  (2)  with  a  suitable  choice  of  f,  the  corrected 
couplings  to  all  higher  partial  waves  are  at  least  three 
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orders  of  magnitude  smaller  than  the  corresponding  uncor¬ 
rected  pss  couplings,  and  the  range  and  magnitude  of  the 
remaining  couplings  are  also  reduced  greatly;  and  (3)  the 
choice  of  the  switching  function  which  produces  this 
cancellation  is  independent  of  continuum  state  character¬ 
istics  and  of  the  type  of  coupling  involved,  but  does 
depend  on  the  initial  discrete  state.  The  results  provide 
strong  evidence  that  the  switching  function  approach  to 
translation  factor  corrections  has  a  sound  physical  basis: 
systematic  reductions  by  several  orders  of  magnitude  are 
not  fortuitous.  This  confirms  earlier  findings,  and 
extends  them  to  additional  discrete  states  and  also  to 
asymmetric  systems.  The  results  are  significant  in  two 
different  contexts. 

2.  Ionisation 

First,  as  was  the  intent  in  reference  2,  the  correc¬ 
ted  couplings  to  continuum  states  can  be  used  to  compute 
physically  sensible  impact  ionisation  cross  sections.  The 
rationale  for  choosing  switching  functions  to  minimise 
continuum  couplings  has  been  described  previously  by 
Thorson  et  al.  c'  a  They  argued  that  (1)  ionisation  cross 
sections  can  be  computed  by  a  first-order  perturbation 
approximation;  (2)  if  the  cross  sections  are  small  (i.e., 
if  coupling  is  weak) ,  one  should  choose  a  zero-order  basis 
set  for  which  the  first-order  couplings  are  as  small  as 
possible;  and  (3)  one  can  look  upon  the  "optimisation" 
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calculations  as  a  kind  of  heuristic  variational  method. 
This  argument  is  intuitive,  and  it  is  not  based  on  any 
formal  variational  principle;  there  is  no  principle  which 
states  that  the  best  zero-order  basis  is  the  one  which 
minimises  the  first-order  T-matrix  (although,  in  the  case 
of  weak  coupling,  this  is  certainly  plausible) .  Note  also 
that  minimising  the  first-order  couplings  does  not  neces¬ 
sarily  minimise  the  first-order  transition  probabilities 
(the  transition  probability  is  essentially  |  /°°  matrix 
element  x  expf-i/^e  dt'/h}  dt|2).  of  course,  if  one 
reduces  a  matrix  element  by  a  factor  of  10 3 ,  it  is  hard  to 
argue  that  the  corresponding  T-matrix  element  is  not 
similarly  reduced.  Finally,  one  may  ask  whether,  in  fact, 
^rst_orc^er  Perturbation  theory  provides  an  adequate  desc¬ 
ription  of  impact  ionisation.  The  whole  question  of  the 
convergence  of  perturbation  expansions  is  extremely  com¬ 
plicated,  but  there  are  reasons  for  believing  that  the 

"off-shell"  contributions  from  the  second-order  terms  may 

57 

well  be  important. 

^1^-hough  the  impact  ionisation  theory  is  based  on 
intuitive  ideas,  the  figures  given  here  and  in  reference  2 
show  that  orders  of  magnitude  differences  will  exist 
between  ionisation  cross  sections  computed  using  uncorrec¬ 
ted  pss  theory ,  and  those  computed  using  the  corrected 

elements.  Tnese  differences  are  exaggerated  by  the 
f^ct  that  the  efficiency  of  ionisation  depends  very 
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strongly  on  the  effective  ionisation  potential  in  the 
region  of  coupling.213'58  m  the  Isa  state  of  H+  for 
example ,  the  ionisation  potential  at  short  distances  is 
significantly  larger  than  the  asymptotic  value,  and  the 
ionisation  cross  section  is  dominated  by  the  contribution 
of  any  couplings  at  distances  greater  than  2  to  4  au.  In 

pss  theory,  these  contributions  are  very  large,  and  almost 
entirely  spurious. 

Figures  25  to  29  suggest  that  ionisation  cross 
sections  for  HeH++  computed  using  the  simple  asymptotic 
switching  functions,  fR  =  ±1,  may  agree  more  closely  with 
those  computed  using  the  "optimum"  f  *s;  indeed,  the  Isa 
and  2sa  cross  sections  should  be  essentially  the  same  in 
each  case.  On  the  other  hand,  substantial  differences 
should  exist  for  "molecular"  states  like  2pa  (ionisation 
of  H ( Is)  by  He++) .  Of  course,  I  think  that  the  more  elab¬ 
orate,  "optimum"  switching  functions  are  the  ones  approp¬ 
riate  for  such  cases. 

3.  Close-coupling 

a.  Comparison  with  pss-type  Theories 

Since  the  examination  of  continuum  couplings  is 
giving  information  about  the  bound  state  f  's,  the  results 
also  have  some  implications  for  the  selection  of  switching 
functions  for  general  close-coupling  problems.  However, 
before  discussing  these,  I  should  first  emphasise  the 
distinction  between  close-coupling  formulations  which 
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correctly  take  account  of  electron  translation  (switching 
function  theories) ,  and  those  which  do  not  do  so  (pss-type 
theories) .  This  distinction  is  important  for  both  formal 
and  practical  reasons. 

I  pointed  out  earlier  that  any  arbitrarily  chosen 
set  of  switching  functions  gives  a  formally  valid  theory 
of  close-coupling.  Moreover,  by  taking  particular  choices 
for  the  f n ' s ,  nearly  all  existing  theories  which  take 
account  of  electron  translation  can  be  accomodated  as 
special  cases  of  such  a  formulation.  Provided  the  basis 
used  is  large  enough,  the  resulting  detailed  cross  sections 
are  9ua*~anteed  to  be  correct  and  independent  of  the  choice 
of  switching  functions.  No  such  guarantees  can  be  given 
for  pss  theory,  because  the  basis  functions  do  not  satisfy 
asymptotic  scattering  boundary  conditions.  Neither  can 
they  be  given  for  modifications  which  take  some  fixed 
point  (e.g.,  the  target,  or  projectile,  nucleus)  as  the 
reference  origin  for  electron  co-ordinates.  This  is  true 
even  though  the  results  may  converge  numerically,  as  the 
number  of  close-coupled  states  considered  is  increased. 

Even  if  we  leave  aside  the  question  of  formal 
validity,  the  uncorrected  theories  are  computationally 
inconvenient.  As  the  figures  show  for  continuum  couplings, 
there  are  many  large,  long-ranged  matrix  elements,  and 
this  is  also  the  case  for  couplings  between  discrete 
states.  The  corrected  couplings  can  be  greatly  reduced  in 
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size  and  range,  even  with  a  very  simple  switching  function. 
This  is  an  important  practical  advantage,  since  computer 
time  is  often  a  limiting  factor  in  close-coupling  studies. 

Recently ,  Winter  and  Lane^  have  described  a  very 
thorough  treatment  of  slow  collisions  in  the  HeH++  system, 
following  an  earlier  study  by  Piacentini  and  Salin.^  In 
both  studies,  the  objective  is  the  calculation  of  the 
total  charge-exchange  cross  section  for  He++  +  H(ls)  col¬ 
lisions.  Both  note  the  problem  of  translation  factors  and 
their  effects,  but  use  the  target  (H)  nucleus  as  the  fixed 
origin  of  electron  co-ordinates.  This  ensures  that  bound¬ 
ary  conditions  are  satisfied  for  all  scattering  states  in 
the  target  channel,  but  not  for  those  in  the  He+  channel. 
The  total  capture  cross  section  can  be  found  from  probab¬ 
ility  conservation .  Winter  and  Lane^  remark  that  they 
did  not  calculate  cross  sections  for  individual  trans¬ 
itions,  because  their  treatment  does  not  fully  take  account 
of  translation  factor  corrections.  They  include  up  to  20 
discrete  states,  and  present  pictures  of  many  of  the 
matrix  elements  obtained  in  this  formulation. 

Choosing  the  proton  as  a  fixed  reference  origin  is 
equivalent  to  taking  fR  =  +1,  for  all  states,  i.e.,  it  is 
not  a  switching  function.  Matrix  elements  coupling  states 
which  are  asymptotically  in  the  He+  channel  behave  like 
pss  matrix  elements;  those  involving  states  in  the  H 
channel  are  asymptotically  well-behaved,  but  are  fairly 
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arge  at  finite  distances.  Using  a  switching  function 
formulation ,  most  of  these  couplings  can  be  greatly 
reduced  in  size  and  range. 

b#  g°mParison  with  other  Switching  Functions 

Finally,  we  may  consider  the  respective  merits  of 
elaborate  "molecular"  switching  functions  (like  those 
determined  here)  ,  and  the  simpler  choices,  f  =  ±1  (to 
which  ours  reduce  as  R~o) ,  for  use  in  close-coupling 
calculations.  Although  the  whole  subject  of  switching 
functions  needs  more  study  before  this  question  can  be 
resolved,  some  relevant  points  are  suggested  by  this  work. 

(1)  As  I  have  pointed  out,  the  behaviour  of  couplings 
to  the  continuum  "selects"  optimum  switching  functions 
which  are  precisely  defined.  Moreover,  they  are  functions 
only  of  the  initial  discrete  state,  and  they  differ  sig¬ 
nificantly  from  the  simple  forms  only  when  the  correspon¬ 
ding  state,  | n> ,  is  genuinely  "molecular."  it  seems  very 
likely  that  the  variation  of  the  switching  finction  is 
reflecting  that  molecular  character. 

(2)  In  a  close-coupling  calculation,  the  basis  is 
truncated  to  a  finite  set  of  molecular  electronic  states, 
and  couplings  to  states  outside  that  manifold  are  ignored. 

In  particular,  couplings  to  the  continuum  are  neglected, 
because  we  know  that  the  ionisation  cross  sections  are 
small.  It  is  certainly  consistent  with  that  neglect  to 
use  a  formulation  in  which  the  theoretical  couplings  to 
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the  continuum  are  made  as  small  as  possible. 

(3)  We  should  recognise  that  it  may  be  possible  to 
reduce  the  corrected  matrix  elements  still  further,  with  a 
more  sophisticated  form  for  f n .  However,  with  the  possible 
exception  of  the  3da  state,  it  is  my  belief  that  it  will 

be  hard  to  improve  much  on  the  spectacular,  three  to  four 
orders  of  magnitude  reductions. 

(4)  It  will  be  interesting  to  compare  the  results 
obtained  by  Winter  and  Lane,20  from  their  close-coupling 
study  of  HeH  ,  with  those  obtained  from  studies  using 
(a)  the  simple  switching  functions,  f  =  ±1,  and  (b)  the 
"optimum"  fn's  reported  here.  In  the  calculations  using 
switching  functions,  it  will  also  be  possible  to  compute, 
and  compare,  the  cross  sections  for  individual  transitions. 


*. 


- 


' 


REFERENCES 


1. 


(a)  W.  R.  Thorson  and  H.  Levy  II,  Phys.  Rev., 

232 . (1969) ; 

(b)  H.  Levy  II  and  W.  R.  Thorson,  Phys.  Rev., 

244  (1969) ; 
ibid. ,  p.  252. 


181, 

181, 


2'  ^  Cm  F*  Lebeda,  W.  R.  Thorson,  and  H.  Levy  II, 

Phys.  Rev.,  A4  ,  900  (1971)  ; 

(b)  V.  SethuRaman,  w.  R.  Thorson,  and  C.  F.  Lebeda, 
Phys.  Rev.,  A8,  1316  (1973). 

3"  175BU958)and  R*  McCarro11'  Proc-  R°y*  Soc.,  A245  , 


4.  Idem,  Adv.  in  Phys.,  11,  39  (1962). 

5"  nq7iilBateS  and  D*  Sprevak'  J*  Phys.,  B3,  14  83 

V  /  v  /  • 


6. 


7. 


8. 

J. 

J. 

9. 

M. 

10. 

A. 

11. 

N. 

12. 

N. 

13. 

J. 

14. 

M. 

S*  3ilS^1969ferman  Russek,  Phys.  Rev.,  181, 

(a)  W.  ^^Tk°^on  end  J*  B*  Delos,  Phys.  Rev.,  A18 , 

(b)  ibid.  ,  p.  135.' 


r  •  JAii  uu.o  *  xr  XI  VS  • 

Rev. ,  A6 ,  709  (1972) ;  y 

B.  Delos  and  W.  R.  Thorson,  Phys.  Rev.,  A6,  720, 

E.  Riley,  Phys.  Rev.,  A7,  626  (1973). 

Messiah,  Quantum  Mechanics,  vol .  II  (Amsterdam* 
North -Ho 11 and  Publ.  Co.,  1965)  p.  747. 

F.  Mott  and  H.  S.  W.  Massey,  Theory  of  Atomic 

i°ns  /  3d  ed.  (Oxford:  Clarendon  Press ,  1965) 

p.  428. 

N.  F.  Mott,  Proc.  Camb.  Phil.  Soc.,  27,  553  (1931). 

P.  Davis  and  W.  R.  Thorson,  Can.  J.  Phys.,  56 
996  (1978) .  — ' 

Born  and  J.  R.  Oppenheimer,  Ann.  der  Phys..  84 
457  (1927) .  y  '  — ' 


172 


K 


* 

. 


•  *  . 


. 


173 


15 


16 

17 

18, 

19. 


20 


21. 


22. 


23 

24 

25 

26, 

27. 


M.  Born  and  K.  Huang,  Dynamical  Theory  of  Crystal 
Lattices,  (Oxford:  Clarendon  Press,  1$ 54 )  ,  ppT 
166-173.  '  ** 

D.  R.  Bates,  H.  S.  W.  Massey,  and  A.  L.  Stewart, 
Proc.  Roy.  Soc.,  A216,  437  (1953). 

A.  Riera  and  A.  Salin,  J.  Phys . ,  B9,  2877  (1976). 

G.  B.  Schmid,  Phys.  Rev.,  A15 ,  1459  (1977). 


R.  D.  Piacentini  and  A.  Salin,  J. 
(1974)  ; 

ibid. ,  B9 ,  563  (1976) ; 
ibid. ,  BIO,  1515  (1977) . 


Phys. ,  B7 ,  1666 


28 


T*  ?^Wl?ter  and  N*  F*  Lane'  Phys.  Rev.,  A17,  66 
(1978)  .  - ' 

C.  F.  Melius  and  W.  A.  Goddard  III,  Chem.  Phys 
Lett. ,  15,  524  (1972) ;  Y 

idem,  Phys.  Rev.  Lett.,  29_,  975  (1972). 

J.  C.  Y.  Chen,  V.  H.  Ponce,  and  K.  M.  Watson,  J. 
Phys. ,  B6 ,  965  (1973) ; 

G.  J.  Hatton,  J.  C.  Y.  Chen,  T.  Ishihara,  and  K.  M. 
Watson,  Phys.  Rev.,  A12,  1281  (1975). 

Y.  Hahn,  Phys.  Rev.,  154 ,  981  (1966). 

A.  V.  Matveyenko,  J.  Phys.,  B9,  1419  (1976). 

M.  H.  Mittleman  and  H.  Tai ,  Phys.  Rev.,  A8 ,  1880 
(1973) .  — 

K.  Taulbjerg,  J.  Vaaben,  and  B.  Fastrup,  Phys.  Rev 

A12,  2325  (1975) .  * 

M.  H.  Mittleman,  Phys.  Rev.,  188,  221  (1969)* 

ibid.,  A9,  704  (1974);  - 

ibid. ,  A10 ,  442  (1974) . 

M.  E.  Riley  and  T.  A.  Green,  Phys.  Rev.,  A4 ,  6' 
(1971)  .  — r 


29.  S.  Cohen,  J.  R.  Hiskes ,  and  R.  J.  Riddell,  Jr.,  Phvs 

Rev. ,  119,  1025  (1960) .  y 

30.  W.  Kolos  and  L.  Wolniewicz,  Rev.  Mod.  Phys.,  35,  473 


' 

■ 


* 


-•.vs 


. 


174 


31.  C.  L.  Beckel, 


B.  D.  Hansen  III,  and  J.  M.  Peek 
.  ,  5_3,  3681  (1970)  . 


J. 


Chem.  Phys 


35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 


F.  T.  Smith,  Phys.  Rev.,  179 ,  111  (1969). 

R'  ff  L*  Kroni9>  Band  Spectra  and  Molecular  Structure, 
(London:  Cambridge  University  Press,  1930) , 


chap.  I. 

34.  T*T  d  mi - 


A.  Messiah,  Quantum  Mechanics,  vol .  I  (Amsterdam: 
North-Holland  Publ.  Co.,  1965)  appendix  B. 

D.  R.  Bates,  K.  Ledsham,  and  A.  L.  Stewart,  Phil. 
Trans.  Roy.  Soc. ,  A246,  215  (1953). 

J.  M.  Peek,  J.  Chem.  Phys.,  £3,  3004  (1965). 

M.  M.  Madsen  and  J.  M.  Peek,  At.  Data,  2,  171  (1971). 

D.  R.  Bates  and  R.  H.  G.  Reed,  Adv.  At.  Mol.  Phys.,  4, 
13  (1968) 

D.  R.  Bates  and  T.  R.  Carson,  Proc.  Roy.  Soc.,  A234, 

207  (1956)  .  '  — — ' 

L.  I.  Ponomarev  and  T.  P..Puzynina,  USSR  Comp.  Math. 
Math.  Phys.,  8(6) ,  94  (1968). 

J.  D.  Power,  Phil.  Trans.  Roy.  Soc.,  A274 ,  663  (1973). 

E.  L.  Ince,  Phil.  Mag.,  6,  547  (1928). 

J.  A.  Stratton,  P.  M.  Morse,  L.  J.  Chu,  and  R.  A. 
Hunter ,  Elliptic  Cylinder  and  Spheroidal  Wave- 
functions-^  (New  York :  JY  Wi  ley  &  sons  ,  19 41). 

K.  Helfrich ,  Theoret.  Chim.  Acta,  21.,  381  (1971). 

J.  H.  Wilkinson,  The  Algebraic  Eigenvalue  Problem, 
(Oxford:  Clarendon  Press,  196  5)  p.  321. - ' 

A.  H.  Stroud  and  D.  Secrest,  Gaussian  Quadrature 
Formulas ,  (Englewood  Cliffs,  N.J.:  Prentice  Hall 


35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 


1966)  . 


■ 


. 


175 


49.  M. 

50 .  An 

51.  D. 

52.  L. 

53.  (a) 
(b) 


A.  Abramowitz  and  I.  A.  Stegun,  Handbook  of 
Mathematical  Functions ,  (Washington,  D.C.:  N.B.S., 
1964)  p.  923. 

article  based  on  the  contents  of  this  section  has 
been  submitted  for  publication  to  J.  Comp.  Phys. 

R.  Bates,  U.  Opik,  and  G.  Poots,  Proc.  Phys.  Soc. 
(Lon. )  ,  66_,  1113  (1953)  . 

183P°l976):eV  and  L*  N*  Somov'  J*  ComP*  Phys.,  20, 


P.  T.  Greenland,  Theoret.  Chim.  Acta,  38,  185 
(1975).  — 

P.  T.  Greenland  and  W.  Greiner,  Theoret.  Chim. 
Acta,  42,  273  (1976) . 


.  S.  D.  Conte  and  C.  de  Boor,  Elementary  Numerical 
Analysis ,  (New  York:  McGraw-Hill,  Inc.,  1972) 
ch  ap .  6  . 


55.  G.  Arfken,  Mathematical  Methods  for  Physicists,  2d 
ed.  (New  York:  Academic  Press,  1970)  chap.  12. 


56. 


J.  Rankin  and  W.  R.  Thorson,  Phys.  Rev. 
(1978)  . 


A00  ,  0000 


57. 


W.  R.  Thorson  and  J.  H.  Choi,  Phys.  Rev.,  A15 . 
(1977)  .  - 


550, 


W.  R.  Thorson,  Phys.  Rev.,  A12,  1365  (1975). 


58. 


. 


- 


;  •  v;  •.  . 


■ 


APPENDIX 


STATE -DEPENDENT  SWITCHING  FUNCTION 
FORMULATION 

A.  Introduction 

The  derivation  of  the  classical  trajectory  equations 
in  chapter  II  assumed  that  the  electron  translation  factor 
can  be  characterised  by  a  single  switching  function  common 
to  all  electronic  states.  It  is  desirable,  however,  to 
extend  the  theory  so  that  each  state  has  its  own,  charac¬ 
teristic  switching  function,  f^.  The  derivation  is  quite 
straightforward,  and  follows  that  of  chapter  II,  section 
C.  It  will  be  shown  that  the  new  coupled  equations  have 
the  same  form  as  before,  but  that  the  definition  of  A(R) 
is  changed.  Although  the  new  correction  matrix  is  no 
longer  Hermitian,  the  equations  still  conserve  probability, 
to  first-order  in  the  nuclear  velocity. 

B.  Definition  of  Basis 
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Fn  =  exp { (im/h)  [v.sn  -  /fc (1-X2) v2dt ' /S] } ,  (A-2a) 


3n  =  %  <fn+*> 


(A-2b) 


and  $n(r;R)  is  a  Born  Oppenheimer  electronic  wavefunction. 

In  chapter  II,  the  states,  *n(r;R),  are  orthogonal, 
but  when  the  switching  functions  are  allowed  to  be  state- 
dependent,  the  overlap  matrix  is  no  longer  diagonal.  Let 
us  define  S_(v)  ,  the  overlap  matrix,  such  that 

Sjn(v)  =  <j|F.*Fn|n> 


=  <j  |exp{  (im/h)  v.  (sn~s  . )  }  |  n> .  (a-3) 

To  first  order  in  v,  we  can  write 

S(v)  =  1  +  v. a ;  S_1(v)  =  1  -  v.o,  (a-4) 


where 


°jn  =  (ira/h)  <j  I  (sn-s  .)  | n>. 


(A-5) 


C.  Classical  Trajectory  Equations 
Assuming  that  the  nuclei  follow  a  classical  path, 
the  wavefunction,  T,  satisfies  the  time-dependent 
Schrddinger  equation. 


he  T  =  ih  3T/3t, 


( A-6 ) 


and  we  can  write  (cf.  equation  (11-42)) 


Proceeding  exactly  as  in  chapter  II,  we  find  that  (cf. 
equation  (II-50a)) 

ih  S(v)  b  =  {h(v)  +  v.  [P(v)+A(v)  ]  }  b,  (A-8) 

where 


hjn(v)  =  <3  I  Fj*Fn  heln><  (A-9a) 

Pjn(v)  =  -ih  <j|Fj*Fn  (VR)j|n>,  (A_9b) 

and 


Ajn(v)  =  (im/h)  <j|F.*Fn  [he,sn]|n>.  (A-9c) 

To  first  order  in  v,  we  can  write 


h(v)  =  £  +  v.n , 


( A-10 ) 


where 


(A-lla) 


njn  =  (im/h)  <j  |  (sn-s.)  he|n>, 


and 


e 


jn 


e  (R)  6. 
n  jn 


(A-llb) 


Thus,  to  first  order  in  v. 


s  1(V) 


h(v)  =  e  + 


v.  (n  - 


a  e} 


(A-12) 


=  £. 

Hence,  if  we  multiply  equation  (A— 8)  on  the  left  by  S~1(v) 
and  keep  only  terms  which  are  first  order  in  v,  we  find 

ih  b  =  {e  +  v.  [P+A]  }  b,  (A-13) 


where 


Pjn  ih  (A- 14 a) 
Ajn  =  (im/h)  <j| [he,sn] |n>.  (A-14b) 

Equations  (A— 13)  are  the  desired  coupled  equations. 
They  have  the  same  form  as  (11-50),  but  the  matrix.  A,  is 
now  non-Hermitian,  and  defined  by  (A-14b) .  I  will  now 
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show  that  probability  is  conserved,  to  first  order  in  the 
nuclear  velocity. 


D.  Conservation  of  Probability 
For  probability  to  be  conserved,  we  require 


d/dt  <T [ T>  =  0. 


(A-15) 


Using  Green's  method,1 


ih  d/dt  <T|T>  =  ih  d/dt  {b+  S (v)  b} 


(A-16 ) 


=  ih  {b+  S(v)  b  +  b+  S(v)  b  +  b +  S  (v)  b}. 


If  we  define 


Pjn  =  (ira/h)  <j |sjn>, 
then,  to  first  order  in  v, 

S  (v)  =  1  +  v.  a  =  1  +  v.  (p  +  p  +  ) 


(A-17) 


(A-18a) 


and 


A 


e  p 


p  £ . 


(A- 18b) 


Using  equation  (A-13) ,  we  can  eliminate  b  and 
from  (A-16) ,  and  obtain 

ih  d/dt  <T|T>  «  b+  {v. (o  e  -  e  a)  +  ih  $.3  (S.S)  + 

R  — 

v.  (A  -  Af)  }  b  ,  (A- 19 ) 

since  both  £  and  P  are  Hermitian  matrices.  But 
A  -  A+  =  £(£  +  £+)  -  (£  +  £^)£ 

=  £  1  "  1  £'*  ( A-20 ) 

and  hence 

ih  d/dt  <T | T>  -  b^  {ih  v. VR(v.a) }  b,  (A-21) 

2 

which  is  of  order  v  . 
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